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I. INTRODUCTION 



A. BACKGROUND 

Surveying is the science and art of measurements which are necessary to 
determine the relative position of points above, on, or beneath the surface of the earth, 
or to establish the points in a specified position. Surveying operations are conducted 
not only on land, but also in the oceans and in space. The measurements of surveying 
consist of distances, horizontal and vertical, and directions. In order to provide a 
framework of survey points whose horizontal and vertical positions are accurately 
known, basic horizontal and vertical control surveys must be performed. A primary 
use of control surveys is for construction of control for a map or chart. The 
fundamental network of points whose horizontal positions have been accurately 
determined is called horizontal control (Schmidt, 1978, p. 122). 

Horizontal control generally is established either by traverse, triangulation, or 
trilateration. Which one is to be used depends on the accuracy required and the factor 
of economy in the selection of survey method. Obviously, there are many degrees of 
precision possible in any measurement because no surveying measurement is exact. 
Each of these methods may be the best one to use for a given purpose. Ordinarily, it is 
a waste of time and money to attain unnecessarily high accuracy. On the other hand, 
if the measurements are not sufficiently precise, faulty survey results are produced. 
Therefore, the best surveyor is not the one who makes the most precise measurements, 
but the one who is able to choose and apply the appropriate measurement with 
precision requisite to the purpose. 

Before 1950 the main framework of a first-order geodetic survey almost always 
consisted of triangulation, which could be replaced by traverse in cases where the 
topography made triangulation impracticable. Today, due to the development of 
electronic distance measuring (EDM) equipment, the first-order control points can be 
established by means of high accuracy traverse [Allan, 1968, p. 370]. Therefore, the 
horizontal control is frequently provided by traverse, especially for surveys in area of 
limited extent, mostly fiat and jungle covered. Traverse in such cases is much more 
economic, convenient, and rapid than other methods and the results are equally 
accurate. 
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In order to achieve high precision of horizontal control points when distributed 
over a large area, first or second-order geodetic surveys are required. These types of 
survey treat the shape of the earth as ellipsoidal and require using the most accurate 
distance and angle measurements. Computation of such a survey is relatively 
complicated, based on long geodetic formulas for computing (with necessary precision) 
the exact horizontal and vertical position of widely distributed points on the earth's 
surface. 

Disregarding ellipsoid shape, a third-order survey is used over earth areas of 
limited extent. In this type of a survey, the earth can be considered as flat and all 
angles are considered to be plane angles. Surveys of this type are used in the 
densification of geodetic control. 

B. OBJECTIVES 

As mentioned above, the traverse method has been used worldwide mostly for 
densification of control stations. However, there are many methods of traverse 
computations. The main objectives of this thesis are to (1) compute a 
closed-connecting traverse and adjust station positions by the Approximation Method 
with the compass (Bowditch) rule and Least Squares Method (adjustment by 
observation equations only), and (2) compare the results of the two methods. 

All computations have been accomplished in the Universal Transverse Mercator 
(UTM) grid coordinates rather than geodetic coordinates. 
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II. TRAVERSES 



A. GENERAL 

1'hc word traverse generally means to pass across. But in surveying, this word 
means the measurement in a specified sequence of the lengths and directions of lines 
between points on the earth whose position may be know or unknown. Traverse is the 
most widely employed method for dcnsification of local horizontal control. Linear 
measurements are made either by direct observation using a tape or EDM equipment, 
or by indirect observation using tachcomctric methods. The angular measurements arc 
made with theodolite or transit. In this thesis, the only traverse operations considered 
arc for angles measured by theodolite and distances measured by precise EDM 
equipment or tapes. 

Two kinds of traverse exist in surveying, an open and a closed traverse. 

I. Open Traverse 

An open traverse normally originates at a point of known position and docs 
not return to the starting point nor docs it terminate on another point of known 
position (Eigurc 2.1). Open traverses should generally not be used because they can 
not be checked for errors. 




Eigurc 2.1 Open Traverse. 



2. Closed Traverse 

Closed traverses can further be sub-divided as, a closed-loop and a 
closed-connecting traverse. 

A closed-loop traverse is one that originates and terminates on a single point 
of known position, thus forming a closed polygon (f igure 2.2). This type of traverse 
provides an internal check on angles but no check on systematic errors in distance. 
Also, if the starting azimuth (between stations 0 and I in Figure 2.2 ) has an error, it 
causes error in orientation of the entire traverse. A closed-loop traverse generally 
should not be used. 




Figure 2.2 Closed-loop Traverse. 



A closed-connecting traverse is one that begins and ends on two different 
points whose horizontal positions have been previously determined by a survey of 
higher or equal accuracy (Figure 2.3). This type of traverse is preferable to all others, 
since computational checks arc possible to detect systematic errors in both distance 
and direction. 

B. ANGULAR AND DIRECTIONAL MEASUREMENTS 

The position of traverse points is determined by the direction and distance from 
the starting point. To obtain the direction by means of azimuth, the horizontal angle, 
or plane angle, must be measured in the field. Also, the determination of vertical 
angles, or zenith distances, may be required to reduce slope distances to horizontal 
distances. 
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Figure 2.3 Closed-connecting Traverse. 

Commonly, both horizontal and vertical angle measurements arc accomplished 
with a transit or theodolite. The theodolite is employed especially for surveys of high 
precision. Two types of theodolite arc a repeating theodolite and direction theodolite. - 
A repeating theodolite reads directly to 20" or T and may be estimated to one-tenth of 
the corresponding direct reading. A direction theodolite usually reads directly to 1" 
and may be estimated to 0.1" | Davis ct ah, 1981, p. 2 1 5 J . In general, a direction 
theodolite is more precise than a repeating theodolite and with it, plane angles arc 
computed by subtracting one direction from another. 

In all types of traverses, the horizontal angles can be measured by one or more 
of the following listed angle measurement methods. 

1. Interior Angle 

Interior angle is the angle measured within a closed figure at the intersection 
of two lines (Figure 2.4). 

2. Deflection Angle 

Deflection angle is the angle measured from the extension of the preceding line 
to the succeeding line (Figure 2.5). Such angles must be identified as right or left to 
express whether the angle is turned to the right or to the left from the prccccding line. 

3. Angle To The Right 

Angle to the right is the clockwise angle measured from the preceding line to 
the succeeding line (Figure 2.6) 
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Figure 2.4 Measuring of Interior Angles. 




Figure 2.5 Measuring of Deflection Angles. 

C. LINEAR MEASUREMENT 

Direct linear measurements may he obtained in traversing by pacing, odometer 
reading, tachcomctry (stadia), subtense bar, taping, and EDM. Of these methods, 
taping and EDM arc most commonly used by surveyors. However, EDM equipment 
has a clear superiority over traditional taping for lines in excess of about 250 meters. 

Distances measured using EDM equipment arc subject to personal, instrumental, 
and natural errors. Personal errors include misreading, improper centering of the 
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Figure 2.6 Measuring of Angles to the Right. 

instrument over the stations, failing to exactly center the null meter, and incorrectly 
measuring meteorological factors and instrument heights. Instrumental errors, 
expressed in terms of the accuracy of the instrument specified by the manufacturer, 
contain two parts. For example, if the accuracy of an instrument is designated as ± 

(10 ppm + 5 mm), the constant error part is + 5 mm, which is independent of the 
distance, and the value of the proportional part is 10 ppm (parts per million) which is a 
function of the distance measured. Constant error is most significant for short 
distances. For very long distances the constant error becomes negligible, but the 
proportional part is important. Natural errors such as refraction arc cause by 
changing of atmospheric conditions along the measured line between the end stations. 

D. ACCURACY 

In survey adjustment, a deviation from the ‘true’ value is considered as an 
observational error and the standard error designates the measure of accuracy of the 
observation. The meaning of an accuracy is then the degree of conformity or closeness 
of a measurement to true value. 

The quality of traverse operations is dependent upon the accuracy of angular and 
linear measurements; thus, in checking the accuracy of traverse two quantities arc 
considered, the angular misciosurc and the linear misclosurc. Although the positional 
closure (relative accuracy) is an indication of the overall quality of the traverse and is 
used for traverse classification, it docs not yield information on the precision of point 
location determined in a traverse [Davis et al., p. 332]. 
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The inherent weakness in a traverse is that the deviation of each measured line is 
determined by a single series of angular observations, further, any error in any angle 
will affect not only the adjoining line but all subsequent lines to a greater or lesser 
extent according to their lengths [Allan, 1968, p. 371], 

Angular misclosure is expressed by standard error of the measured angle times 
the square root of the number of measured angles. 

Linear misclosure is commonly expressed as a ratio of total misclosure to total 
length of traverse. 

Finally, some of the most significant features of traverse classification by the U.S. 
Federal Geodetic Control Committee (1984) are shown in Table I. 

E. ADJUSTMENTS 

Adjustment of a traverse is carried out to ensure consistency within the known 
positions of the originating and terminating stations and to remove inconsistencies in 
observed angles and distances to compensate for random errors. For a more precise 
extended traverse, adjustments made on the basis of least squares arc preferred. But a 
traverse of limited extent can be adjusted by simple approximation methods. 

1. Approximation Methods 

There are four methods for traverse adjustment by approximation. 

a. Arbitary Method 

This method does not conform to a fixed rule. Rather, the linear error of 
closure is distributed arbitarily according to arbitary preference of the surveyor. 

b. Transit Rule 

Transit rule is better for adjustment of the traverse where the angles are 
measured with greater accuracy than distances, and is valid only when the traverse lines 
are parallel with the grid system used for the traverse computations. Corrections are 
made by the following rules: the correction in latitude for any station is equal to the 
multiple of latitude in that section and total closure in latitude divided by the sum of 
all latitudes in traverse, and the correction in departure is equal to the multiple of 
departure in that section and total closure in departure divided by the sum of all 
departures in the traverse [Davis et al., 1981, p. 323]. 

c. Compass or Bowditch Rule 

This method is suitable for adjustment of the traverse where the angles and 
distances are measured with equal precision and uses the following rules: the correction 
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TABLE I 

TRAVERSE CLASSIFICATION 



Order 

Class 


First 


Second 

I 


Second 

II 


Third 

I 


Third 

II 


Station spacing not 
less than (km) 


10 


4 


2 


0. 5 


0. 5 


Minimum number of 
network control 
points 


4 


3 


2 


2 


2 


Theodolite least 
count 


0 . 2 " 


1 . 0 " 


1 . 0 " 


1 . 0 " 


1 . 0 " 


Direction number of 
positions 


16 


8 or 12 


6 or 8 


4 


2 


Standard deviation of 
mean not to exceed 


0. 4" 


0. 5" 


0 . 8 " 


1 . 2 " 


2 . 0 " 


Rejection limit from 
the mean 


4 » 


5" 


5" 


5" 


5" 


Azimuth closure at 
azimuth check point 
( second of arc) 


1 . 7VN 30VN 


4. 5VN 10. OVN 


12 VN 



Position closure 0. 04,/K 0. 08 Vk 0. 2>/k 0. 4,/K 0. 8 Vk 
after azimuth or or or or or 

adjustment* 1: 100000 1: 50000 1: 20000 1: 10000 1: 5000 

N = route distance in km 
K = number of segments 

*The expression containing the square root is designed 
for longer lines where higher proportional accuracy is 
required and the results are in meter. 

The closure ( e. g. , 1:50,000) is relative error of 
closure. 



in latitude for any station is equal to the multiple of the length in that section and total 
closure in latitude divided by the total length in the traverse, and the correction in 
departure is equal to the multiple of the length in that section and total closure in 
departure divided by the the total length in the traverse [Schmidt, 1978, p. 150]. 
d. Crandall Method 

Crandall method is a rather complicated procedure which is more rigorous 
than either the compass or transit rules but suitable for adjusting traverses where the 
linear measurements contain larger random errors than the angular measurement. In 
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this method, the angular error of closure is first distributed in equal portions to all of 
the measured angles, then linear measurements are adjusted by using a weighted least 
squares procedure [Brinker, 1977, p. 22S], 

2. Adjustment by Least Squares Method 

The method of least squares adjustment is based upon the theory of 
probability; it simultaneously adjusts the angular and linear measurements to make the 
sum of the square of the residuals (error) a minimum [Brinker, 1977, p. 228], This 
method can be used for any type of traverse. Because of the availibility of fast 
computing devices at the present time, the least squares method is being widely used. 
Further, the least squares solution has the advantage that it determines, quite 
objectively, a unique solution for a given adjustment problem [Clark, 1973, p. 121]. 

In general, adjustment is needed whenever there are redundant observations 
(more observations than are necessary to solve the required unknowns). As an 
example, to determine the angles of a plane triangle, only two observed angles are 
required because the third angle can be obtained by subtraction from 180°. When 
three angles are observed, the sum of them will not be equal 180° due to error in 
measurements. Therefore, these three angles should be adjusted to fit the functional 
model. 

The redundancy may be interpreted to mean that among n observations there 
exist r conditions or functions (n > r) that must satisfy the model. 

Let n be a number of observations and nQ the minimum number of 
observations to find the uniquely solution in the model, then redundancy or degree of 
freedom in the statistic, r, is 

r = n - n 0 (2.1) 

Consequently, there are r redundant observations, which can also give a 
solution. To detect the error in each observation, the best estimated or the most 
probable value must be defined because the true value is not known exactly. 
Statistically, the best estimated value of a group of repeated observation is the average 
(arithmetic mean). 

Once the difference between observed value (X Q ) and estimated value (X £ ) is 
determined, the adjusted value (X a ) is obtained through a least squares solution, then 
the residual (v) can be expressed as 



17 



V 



( 2 . 2 ) 




and 



X a = X e + dx 



(2.3) 



where dx is the correction to estimated value to obtain the adjusted value. 

The least squares adjustment method is based upon the criterion of the sum of 
the squares of the observational residuals must be minimum. 

When observations are considered as uncorrelatcd and of equal precision (with 
identity weight matrix), the least squares condition can be expressed as 

<I> = + = S v i^ = minimum (2.4) 

or in matrix form 

(j> = yTy = minimum (2.5) 

where V is the vector of residuals. 

In uncorrelatcd observations with unequal precision, such as distances and 
angles [Mikhail, 1981, p. 68], the Equations 2.4 and 2.5 become 

0 = WjVj 2 + W 2 V 2 2 + ... + w n v n 2 = £ w i v i 2 = minimum (2.6) 

or in matrix form 

( I> = V^WV = minimum (2.7) 

where wj is the i 1 * 1 clement of the diagonal weight matrix W and vj is the residual 
associated with the corresponding i 1 ^ observation. 

Generally, the relative weights are inversely proportional to variance, thus the 
weight matrix is the inverse of cofactor matrix, Q (when it is square and nonsingular) 
and defined as 



W = Q' 1 



( 2 . 8 ) 



18 



where the elements of cofactor matrix Q are 

% = (2-9) 

and 

qjj = CTij/o-Q 2 (2.10) 

where <Tj~ the variance of the i tF observation, <7- is the covariance between the i tF and 
j in observations, and <7 q Z is variance of unit weight [Mikhail, 1981, p. 67]. 

For the case of uncorrelated weight observations, the cofactor matrice will be 
diagonal with all off-diagonal elements being equal to zero, thus the diagonal elements 
of weight matrix in this case are 

w ii = l '% = V^ii 2 (2- 11 ) 

Generally, there are two types of equation in least squares adjustment: 
condition and observation equations. Condition equations include one or more 
observations but observation equations include parameters and only one observation. 

The condition as well as the observation equations involved in an adjustment 
problem can be linear or nonlinear. However, least squares treatments are generally 
performed with linear functions, since it is rather difficult and often impractical to solve 
nonlinear models [Mikhail, 1978, p. 108]. Consequently, whenever the equations in the 
model are originally nonlinear, they have to be linearized. A method of series 
expansion, especially Taylor's series, is often used to obtain linear equations. Only the 
zero and first-order terms are used and all other higher-order terms are neglected. 

Thus, the linearized form for the general case of m functions of n variables becomes 

F = F ° + J mn Ax (2-12) 

where is the zero-order terms, when x = xP, and J nm is a Jacobian matrix of 
coefficients of first order of n variable (Appendix A). 

The choice between the observation equations (indirect observation) and 
condition equations (observation only) techniques will depend mainly on the 
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mathematical model of the problem to be solved. However, the final answers are 
always the same. 

a. Adjustment by Observation Equations 

The method of adjustment by observation equations is performed with both 
observations and parameters. Therefore, the number of equations is equal to the 
number of observations. Using the example at the beginning of this section, if three 
measured angles and their residuals in a plane triangle are assummed to be ct, P, y, Vj, 
V 2 , and v 3 , respectively, and the adjusted value of those angles is x a j , x a 2 > an< ^ x a 3 
then the three condition equations in the normal form (zero at the right-hand side) are 
(a + v i)-x a i = 0 

(P + v 2> * x a 2 = 0 
(Y + v 3 )-x a3 = 0 



by using Equation 2.3 for adjusted values 

V 1 ' dx l = x eU a 
v 2 ' dx 2 = x e 2 ' P 
v 3 • dx 3 = x e3* y 

letting 





V 1 




-1 


0 


0 




dx-. 


V = 


v 2 

v 3 


/ B = 


0 

0 


-1 

0 


0 

-1 


, X = 


dxi 

dx§ 



then, expressed in matrix notation as 



x 



x 



el 

e2 

e3 






Y 



V 31 + B 33 X 3 j = F 31 and the general form of the adjustment of observation 
equations become 



V nl + ®nu^ul 



= F 



nl 



(2.13) 



where V is an n by 1 vector of residuals; B is an n by u matrix of coefficients; X is an u 
by 1 unknown vector in which u is the number of parameters; and f is an n by 1 vector 
and equals P - X Q , in which P is evaluated using the estimated value X £ . 
b. Adjustment by Condition Equations 

The method of adjustment of condition equations only has no parameters 
included in the condition equations. Thus, the number of condition equations is equal 
to the number of redundancies. From the example above, n= 3, nQ = 2, and r= 1, 
which is the number of condition equations being set. In this case, because the sum of 
interior angles must equal 180°, the single condition equation can be expressed as 
a + P + Y + v j + V 2 + v 3 = 180° 
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v 1 + v 2 + v 3 = 1 80° - a - p - y 



A = | 1,1,1 | , 



V = 



v ] 

v- 



Then, the general form of this technique is 



| 180-a-p-y | 



A rn V nr = 



rl 



(2.14) 



When the conditions are originally linear, the vector F is usually written in terms of the 
given observations as 



rl 



- P d - A rn X o nl 



(2.15) 



where A is the coefficient matrix V., P is a constant term (see Section Il.E.a), X Q is 
observed values, r is redundancy, and n is a number of observation [Mikhail, 1976, p. 
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III. DATA ACQUISITION AND REDUCTION 



A. DATA ACQUISITION 

Taverse data used in this thesis were obtained from field work accomplished from 
25 September thru 9 October 1972 by CAPT Glen R. Schaefer, NOAA Corps, and Mr. 
Jim D. Shea, National Ocean Service (NOS), utilizing traverse methods in Pinellas 
County, Florida. Only the first 15 of 40 occupied stations and three intersection points 
will be used for analysis (Figure 3.1). The two pairs of known stations for this 
closed-connecting traverse are shown in Table II. 

The known stations were observed by the US Coast and Geodetic Survey (now 
NOS) and adjusted by the National Geodetic Survey (NGS). Station Turtle 2 is of 
first-order and the other three stations are of third-order. 

The horizontal angles measured by the method of angles to the right (Table III) 
were turned with a Wild T-2 theodolite, according to the specifications of third-order 
class I traverse, by starting at station Tomlinson and using Egmont Key Lt. House for 
a backsight. The traverse was closed on Turtle 2 with a check azimuth to Madeira 
Beach Tank. 

The slope distances were measured in the field with a Model 76 Geodimeter in 
feet and corrected for temperature and pressure. Distance measurement by the Model 
76 Geodimeter are reported to have an accuracy in the temperature range of -20°C to 
+ 50°C of ±(1 ppm + 1 cm) with a resolution of 1 mm. [Schmidt, 1978, p. 116]. all 
observed distances were converted to meters and reduced to horizontal by the 
procedure given later in this chapter. Finally, geodetic distances are reduced to grid 
distances by applying the scale factor correction (Table IV). 

B. DATA REDUCTION 

1. Computation of Ellipsoidal Distances 

For the requirement of high precision in the first-order traverse, the measured 
distances obtained by EDM equipment are first corrected for atmospheric conditions 
and then reduced to the ellipsoid (Figure 3.2) by the equations 

S = 2 R a Sin' 1 (S 0 /2R a ) (3.1) 
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A = known station 
O = traverse station 
IP = intersection point 



Nor ro Scale 




Figure 3.1 Traverse Layout. 
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TABLE II 








DATA OF KNOWN STATIONS 




Station 


Station 


Grid coordinate ( m. ) 


No. 


Name 


X/Easting 


Y/Northing 


0 


Egmont Key Lt. House 


326214. 833 


3054014. 779 


1 


Tomlinson 


329400. 420 


3063485. 483 


15 


Turtle 2 


325348. 292 


3076179. 297 


16 


Madeira Beach Tank 


322698. 706 


3076362. 405 



s 0 = { [m 2 -(h 2 -h,)] / (1 + hj/RgXl + h 2 /R a ) } ^ 2 



(3.2) 



1/R a = {(Cos 2 a)/M} + ((Sin 2 a)/N} (3.3) 



M = a(l-e 2 )/(l-e 2 Sin 2 <p) 3 / 2 (3.4) 



N = a / (1 - e 2 Sin 2 <p) 1 / 2 (3.5) 

where m is the measured distance, S is the ellipsoidal distance, Sq is chord distance, hj 
and h 2 are the heights above the ellipsoid at each end of the line, R a is the radius of 
curvature of the chord distance m at an azimuth a, a is geodetic azimuth clockwise 
from north, M is radius of curvature in the plane of the meridian, N is radius of 
curvature in the prime vertical, a is a semimajor axis of the reference ellipsoid, e is its 
eccentricity, and <p is the mean latitude of that line. [Torge, 1980, pp. 48,49,50,51,180] 
2. Computation of Geodetic Distances 

For a lower-order traverse, the measured distance can be reduced to mean sea 
level (MSL) or geoid only. Because the error of only 1 ppm in length results from an 
error of 6 m in separation between MSL and the spheroid [Bomford, 1980, pp. 42, 342, 
345]. 
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TABLE III 



OBSERVED HORIZONTAL ANGLES 



A) For Traverse 



At Stn. Stn. 


Angles 


No. 


Name 


D 


M 


S 


0 


Egmont Key Lt. House 








1 


Tomlinson 


105 


21 


25 


2 


TN-01 


243 


39 


18 


3 


TN-02 


168 


10 


15 


4 


TN-03 


59 


55 


56 


5 


Ruscue 


291 


28 


29 


6 


RE-01 


160 


43 


42 


7 


RE-02 


269 


55 


50 


8 


RE-03 


92 


43 


29 


9 


RE-04 


178 


22 


51 


10 


RE-05 


182 


31 


19 


11 


RE-06 


196 


54 


42 


12 


RE-08 


168 


51 


46 


13 


RE-09 


161 


06 


42 


14 


RE-10 


236 


58 


14 


15 


Turtle 2 


78 


37 


24 


16 


Madeira Beach Tank 









<7 

( ±) 



B) For Intersection Points ( IP) 



IP #1 * 





At Traverse 


Angles 


Back Stn. 


Stn. 


D 


M 


S 


RE-01 


RE-02 


196 


23 


46 


RE-04 


RE-05 


184 


34 


22 


RE-05 


RE-06 


186 


29 


15 


RE-06 


RE-08 


34 


43 


31 


RE-08 


RE-09 


2 


41 


22 


* * 


RE-03 


RE-04 


175 


35 


01 


RE-04 


RE-05 


97 


12 


06 


RE-05 


RE-06 


2 


32 


22 


k k k 


RE-08 


RE-09 


194 


01 


22 


RE-09 


RE-10 


232 


11 


34 



* 

•k k 



St. Petersburg BCH CO Tank 
St. Petersburg BCH St. Johns CH Tower 
*** Bay Pines Veterans Administration Hosp. 
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TABU: IV 

MFASURFD DISTANCES AND GRID DIS 


I'ANCFS 




At St n. 


Measured 


Distances 


Grid Distances 


a 


To Stn. 


ft 


m 


ft 


m 


( ±m) 


1 


2300. 98 


701. 340 


2300. 91 


701. 318 


0. 017 


2 


2605. 41 


794. 131 


2605. 28 


794. 090 


0. 018 


3 


4452. 48 


1357. 119 


4452. 31 


1357. 068 


0. 024 


4 


709. 56 


216. 274 


709. 54 


216. 267 


0. 012 


5 


5050. 37 


1539. 356 


5050. 21 


1539. 307 


0. 025 


6 


6620. 54 


2017. 945 


6620. 33 


2017. 880 


0. 030 


7 


1655. 17 


504. 497 


1655. 12 


504. 481 


0. 015 


8 


1605. 85 


489. 464 


1605. 80 


489. 448 


0. 015 


9 


2114. 68 


644. 556 


2114. 61 


644. 535 


0. 016 


10 


2362. 13 


719. 979 


2362. 05 


719. 956 


0. 017 


11 


1360. 16 


414. 578 


1360. 10 


414. 561 


0. 014 


12 


5060. 98 


1542. 590 


5060. 41 


1542. 415 


0. 025 


13 


6172. 26 


1881. 309 


6171. 67 


1881. 128 


0. 029 


14 

15 


6664. 72 


2031. 411 


6664. 51 


2031. 346 


0. 030 



m 


= 


measured distance 






So 


= 


chord distance 




S 


S 


= 


arc distance on 


\ 








ellipsoid 




So y 


h , 


= 


elevation above 










ellipsoid at 










station 1 




\ A* 


'h 




elevation above 










ellipsoid at station 


2 






= 


radius earth curvature 








along measured line 




0 



Figure 3.2 1 Ellipsoidal Distance. 

However, the process of reduction requires three steps: (1) correct slope 
distances to horizontal distances, (2) reduce horizontal distance to geodetic distances, 
and (3) change the geodetic distances to grid distances. 

1 he slope distance data used included correction for atmospheric conditions. 
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In a plane survey, such as a traverse, there are two considerations for the 
reduction of slopojdistances to horizontal distances, a short slope distance and a long 
slope distance. 

a. Slope Reduction for Short Distances 

Short slope distances ( < 2 ini or 3.3 km) measured by using IiDM 
instruments separate from the theodolite, can be reduced to horizontal with a simple 
trigonometry process as 

D = m cos 0 (3.6) 

where D is the horizontal distance, m is the slope distance, and 0 is the vertical angle 
(Figure 3.3). 




Figure 3.3 Slope Reduction for Typical Triangle. 

The horizontal distance also may be determined by using the difference in 
elevation between the two ends of the line. The horizontal distance is 

D = (n^-d 2 ) 1 ' 2 (3.7) 



in which d is the difference in elevation between the two end points. The heights of the 
EDM instrument and reflector above the survey mark must be observed, and d 
becomes 
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d = (E 1 + i 1 )-(E 2 + o 2 ) 



(3.8) 



where Ej and E? are the elevation at each end of the line respectively, I j is the height 
of the instrument, and 0 2 is the height of the reflector. Then, expanding the right side 
of equation 3.7 with the binomial theorem yields 

D = m-(d 2 /2S + d 4 /8S 3 + ...) (3.9) 

The quantity in the parenthesis is designated slope correction. For 
moderate slope the first term is usually adequate. When the slope distances and 
vertical angles are obtained by separated EDM equipment from theodolite, the 
correction of the vertical angle must be determined. The corrected vertical angle 0y is 

e T = e 0+ A0 (3.10) 

where 0 q is an observed vertical angle by theodolite, and A0 is 

A0 = H Cos 0 Q / S Sin 1" (3.11) 

and 



H = (H r -H t )-(H r H e ) (3.12) 

where H r is the height of the reflector, H t is the height of the target, Hj is the height of 
the EDM, and H e is the height of the theodolite [Davis et al., 1981, pp. 103-104]. 
Equations 3.11 and 3.12 are not needed when the slope distances and vertical angles 
are obtained simultaneously by using an EDM transmitter built into a theodolite, 
b. Slope Reduction for Long Distances 

Slope reduction for long distances (> 2 mi or 3.3 km) involves using 
vertical angles affected by curvature and refraction. By assuming a mean radius for the 
earth of 3959 mi or 6371 km, then the curvature correction (C), expressed as an angle 
in seconds, is 4.935" per 1000 ft or 16.19" per km and the horizontal distance, D, is 

D = m Sin(9O°-0-C) / Sin(90° + C) (3.13) 
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for 



6 = (Y + P)/2 (3.14) 

where y and P are the vertical angles at each end of the measured line [Davis et al., pp. 
106-107], 

When a single vertical angle (y) is observed, 0 is the corrected vertical angle 
for combined results of curvature and refraction (C&R), then, 0 is y + (C&R)'\ The 
C&R correction is 4.244" per 1000 ft or 13.925" per 1000 m. The correction of C&R 
will be positive when the vertical angle is an elevation angle and negative in the case of 
a depression angle [Davis et al, 1981, p, 108]. 

3. Reduction of Horizontal Distances to Geodetic Distances 

The horizontal distance at same elevation above the geoid, must be reduced to 
a geodetic distance. This can be done by the equation 

D' = (R)(D)/(R+E) (3.15) 

where D' is the geodetic distance, R is the mean radius of the earth's surface at that 
section, D is the horizontal distance at elevation E above the geoid [Davis et al., 1981, 
p. 107], 

C. GRID DISTANCES 

The traverse computation, based on the UTM grid coordinate system, requires 
the reduction of geodetic distances to the plane of the projection by applying the 
projection scale factor and grid scale constant. Scale factor can be obtained from a 
graph or from a rigorous formula [Department of the Army Technical Manual, 1958, 
pp. 3,4,9,17]; i.e., 

k = k Q [ 1 + (XVIII) q 2 + 0.00003 q 4 ] (3.16) 

where k is the scale factor at scale working on the projection, kQ is the central scale 
factor which is an arbitary reduction applied to all geodetic lengths to reduce the 
maximum scale distortion of the projection (for UTM, kQ = 0.9996), and values for q 
and (XVIII) are obtained by the formulas which shown in Table V. 
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TABLE V 






SPECIFICATION OF PARAMETERS 


XVIII 


= 


1 + e' Cos 2 <p . 1 . 10 12 

2 v 2 






a 


V 




( 1 - e 2 Sin 2 (p ) 1//2 


e 2 


= 


( ecentricity ) 2 = ( a 2 - b 2 ) / a 2 


e 2 ' 


= 


e 2 / (1 - e 2 ) 


q 


= 


0. 000001 E' 


E 


— 


grid easting = E' + 500,000 when point is 
east of central meridian, 500,000 - E 
when point is west of central meridian 


V 


= 


radius of curvature in the prime vertical 


a 


= 


semi-major axis of the ellipsoid 


b 


— 


semi-minor axis of the ellipsoid 


<P 


— 


latitude 
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IV. TRAVERSE COMPUTATION AND ADJUSTMENT 



Linear measurements and angles must be checked by computation to determine 
the position of traverse stations and whether the traverse meets required precision. 
Traverse station coordinates are usually expressed in terms of geographical coordinates 
or rectangular coordinates such as those based on an UTM projection. Traverse 
computations are usually done in rectangular coordinates because of the ease of 
computation. In this thesis, only closed-connecting traverse computations in UTM 
coordinates were used. When specifications were satisfied, the traverse was adjusted 
for perfect geometric consistency among angles and lengths. 

A. DATA PROCESSING 

1. Set up of data base 

Two files on IBM 370/3033AP main frame computer at NPS were established. 

2. Modification of an Existing Program 

The Fortran program TRAVADJ, originally written by LCDR Saman 
Aumchantr, RTN (1984) in Watfiv language for computing and adjusting the traverse 
station coordinates, was modified to be able to handle the distances reduction 
processes. 

3. Writing a New Program 

A Fortran program INDTRA was written for computing and simultaneously 
adjusting traverse station and intersection point coordinates by least squares 
observation equations method. 

B. COMPUTATION OF STARTING AND CLOSING AZIMUTHS 

The directions of lines by means of azimuth are used for traverse computation, 
because sines and cosines of azimuth angles automatically provide correct algebraic 
signs for latitudes and departures. 

The terms latitude and departure are widely used in rectangular coordinate 
calculations of surveying. The latitude of a line is its projection on the reference 
meridian, which differs from geographic latitude. The departure of a line is its 
projection on the east-west line perpendicular to the reference meridian. In traverse 
calculations, north latitudes and east departures are considered plus; south latitudes 
and west departure, minus. 
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Latitudes are also sometimes termed ‘northings’ and ‘Y dilFerences’ (AY); 
departures are similarly called ‘eastings’ and ‘X differences’ (AX). 

Because the closure angle of traverse can be checked by azimuth of each 
consecutive line, starting and closing azimuths have to be determined in the first step 
of computation and azimuths can be computed from a pair of known coordinate 
station positions at the two ends of the traverse (Figure 3.1) 

The azimuth of the line from A to B, ct A g, measured clockwise from north, is 
determined by the equation 

a AB = Tan-*(AX/ AY) (4.1) 



for 



AX - X„. X A 



(4.2) 



and 



AY - Y B- y a < 4 ' 3 > 

where X A and Xg are the grid easting coordinates, and Y A and Yg are the grid 
northing coordinates of stations A and B, respectively (Figure 4.1). The quadrant of 
the azimuth of line AB, a A g, is dependent upon the sign of AX and AY (Table VI). 
The back azimuth dg A (the azimuth from B to A) is obtained by adding 180° to the 
forward azimuth « A g. 

The length of the line AB (denoted as S A g or S) can be determined by the 
Pythagorean theorem or by one of the trigonometric relationships 

S = AX / Sin a (4.4) 



or 



S = AY / Cos a 



(4.5) 
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Figure 4.1 Azimuth Computation. 



TABLE VI 

QUADRANT OF AZ1MUT1 1 



Quadrant of azimuth measured Sign of Sign of 
clockwise from north AX Af 



0 to 90 + + 

90 to 180 + 

180 to 270 

270 to 360 - + 



Substituting data from Table 11 into liquations 4.2 and 4.3, the AX and AY 
between Fgmont Key Lt. House and Tomlinson arc computed as AX = 3185.585 m 
and AY = 9470.704 nr. The azimuth of' the line from Fgmont Key Lt. House to 
Tomlinson is 18° 35' 27.6" (Equation 4.1). Similarly, the azimuth Trom Turtle 2 to 
Madeira Beach lank is computed to be 273° 57' 12.0". These starting and closing 
azimuths will be used Tor computing the coordinates of each traverse station and lor 
checking the angular error. 
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C. COMPUTATION OF TRAVERSE STATION COORDINATES 

Computation, of traverse station coordinates is the reverse process of finding 
azimuth and distance from coordinates. Therefore, the rectangular coordinates for 
each closed-connecting traverse station can not be computed unless forward azimuth 
and distance from the previous station are known. 

The azimuth is reckoned clockwise from north and obtained by 

a jk = a- + 180° + Pj (4.6) 

where oij k is the forward azimuth from station j to station k, a- is the forward azimuth 
from station i to station j, and Pj is the horizontal angle at station j for j values of 1 to 
n, where the previous station i = j - 1 and the next station k = j + 1. The number j 
will increase from 1 (which designates the starting known station of the traverse) to 
number n (which was the last occupied and known station of the traverse). 

Departures and latitudes are then computed by using Equations 4.4 and 4.5 
which are rewritten as 

(4.7) 

and 



^^jk Sjk^ n a jk 



AY jk - S jk Cos « jk (4.8) 

where Sj k is the distance between stations j and k. The coordinates of all other 
traverse stations can be determined by adding successive departures (AXj k ) and 
latitudes (AYj k ) to the previous station's X and Y coordinates, respectively. 

Using the data in Tables II, III, and IV, the azimuth and coordinate 
computations at the first station are shown here by using Equations 4.6, 4.7, and 4.8. 
The result for other stations can be seen in Table VII. 

Computing the azimuth from station 1 to 2 

Az = 18°35'27.6" + ISO 0 + 105°21'25" = 303°56'52.6" 

Computing the coordinates at station 2 

X-easting = 329400.420 + [701.318 Sin(303°56'52.6")] 

= 328818.645 m 
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Y-northing = 3063485.483 + [701.3 18 Cos(303°56'52.6")] 

= 3063877.126 m 

When the coordinates of all stations have been computed, they are still 
unadjusted coordinates and can then be adjusted by one of the two methods of Section 
II.E (approximation and/or least squares method). 



TABLE VII 

UNADJUSTED TRAVERSE STATION POSITIONS 



Stn. Angles Forward Dist. Unadjusted 

Azimuths (m) Coordinates 

D M S D M S X(m) Y(m) 











18 


35 


28 






1 


105 


21 


25 








* 329400. 420 


* 3063485 . 483 










303 


56 


53 


701 . 318 




2 


243 


39 


18 








328818. 645 


3063877 . 126 










7 


36 


11 


794 . 090 




3 


168 


10 


15 








328923. 709 


3064664 . 236 










355 


46 


26 


1357 . 068 




4 


59 


55 


56 








328823. 700 


3066017. 613 










235 


42 


22 


216 . 267 




5 


291 


28 


29 








328645 . 029 


3065895 . 760 










347 


10 


51 


1539. 307 




6 


160 


43 


42 








328303. 493 


3067396 . 699 










327 


54 


33 


2017 . 880 




7 


269 


55 


50 








327231 . 464 


3069106 . 259 










57 


50 


23 


504 . 481 




8 


92 


43 


29 








327658. 538 


3069374 . 789 










330 


33 


52 


489. 448 




9 


178 


22 


51 








327418. 001 


3069801 . 054 










328 


56 


43 


644. 535 




10 


182 


31 


19 








327085. 512 


3070353.210 










331 


28 


2 


719. 956 




11 


196 


54 


42 








326741. 616 


3070985 . 723 










348 


22 


44 


414. 561 




12 


168 


51 


46 








326658. 106 


3071391 . 785 










337 


14 


30 


1542. 415 




13 


161 


6 


42 








326061 . 428 


3072814 . 113 










318 


21 


12 


1881 . 128 




14 


236 


58 


14 








324811 . 349 


3074219 . 797 










15 


19 


26 


2031 . 346 




15 


78 


37 


24 


273 


56 


50 


325348. 180 


3076178 . 924 



* Coordinates for station 1 were known and held fixed. 
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D. ADJUSTMENT BY APPROXIMATION METHOD 

In this thesisT the method of Compass or Bowditch rule was used to adjusted the 
data in Tables III and IV. Thus, the first step is to determine the angular error of 
closure and adjust the angles to obtain the proper closing azimuth (closed azimuth at 
fixed stations). 

1. Angular Errors of Closure 

In the closed-connecting traverse, when there are n stations of observed 
horizontal angles, n-1 lines will be measured. An angular error in traverse can be 
checked and obtained at the last station by comparing the computed azimuth and 
closing azimuth at the known station. 

The closing azimuth computed (from the known station coordinates at the 
traverse end) at the station 15 is 273° 57' 12.0". But because of error in measurement, 
the azimuth computed through the traverse at this station is 273° 56' 49.6", which is a 
difficiency of 22.4". This amount of angular error in 15 observed stations meets the 
limit for allowable error for a third-order class I traverse (allowable error from Table I 
is 38.7"). 

The average correction (Table VIII) is distributed uniformly over all the 15 
traverse angles (Table III). 

2. Linear Errors of Closure 

When all angles have been corrected, the process of calculating the improved 
coordinates of all traverse stations may be done. The check on angular closure for a 
closed traverse does not guarantee that the entire survey is correct, because there can 
be considerable errors in the linear measurement of individual lines. Such errors may 
not show up in the angular check. In order to check the closure of the traverse, it is 
necessary to determine linear error. 

The linear error (LE), the departure error (5x), and latitude error (Sy) in a 
traverse are determined by equations 

LE = [ (Sx)2 + (5y)2 ]E 2 (4.9) 



5x = 



GE n - GE n 



(4.10) 
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(4.11) 



6y - GN n ■ GN n 



where 5x and Sy are the traverse closure in departure and latitude, GE n and GE n ' are 
the known and computed grid easting, and GN n and GN n ' are the known and 
computed grid northing at the closing station, respectively. By substituting the fixed 
and computed values of grid easting and northing in Equations 4.10 and 4.11 for the 
data of Tables II and VII 

6x = 325348.292 - 325348.180 = + 0.112 m 

Sy = 3076179.297 - 3076178.924 = + 0.373 m 
LE = [(0. 1 1 2) 2 + (0.373) 2 ] 1/2 = 0.389 m 

The relative error of closure provides a better assessment of the quality of a 
traverse than the linear error of closure. Therefore, it is common practice to calculate 
the relative error of closure, which is the linear error of closure divided by total 
distances of traverse, and to express the result in the form of a ratio with unity as the 
numerator. For the data of Tables II and VII, this computation is 0.389 / 14853.800 
or 1 : 38,185. 

Using the Compass or Bowditch rule, the computed traverse closures and 
corrections were obtained (Tables VIII and IX) and then the adjusted station 
coordinates (Table X) were computed. 



E. LEAST SQUARES ADJUSTMENT BY OBSERVATION EQUATIONS 

The adjustment by observation equations, shown in general form by Equation 
2.13, can be done more directly than the adjustment by condition equations. To 
achieve this, the explicit expression for the residual vector V from Equation 2.13 is 
substituted in Equation 2.7 to obtain the following equation: 

O = (F - BX) T W (F - BX) (4.12) 

= (F T - B T X T )(WF - WBX) 

= X T B T WBX - X T B T WF + F T \VF - F T \VBX 
= X T B T WBX - 2F T \VBX + F T WF 

where X^B^WF = F^WBX are scalar quantities. 
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TABLE VIII 
TRAVERSE CLOSURE 



I) Angular error 



Known azimuth at last station 

Computed azimuth 
Angular error 

Angular correction per station 



273° 57' 11.95" 

273° 56' 49.58" 
22. 37" 



II) Linear error before adjusting azimuths 



Known Coordinate 
at Turtle 2 



X(m) 

325348. 292 



Y(m) 

3076179. 297 



Computed Coordinate 

at Turtle 2 325348.180 3076178.924 

Error 5x = + 0. 112 Sy = + 0. 373 



Linear error of closure 
Total distances 
Relative error of closure 



0. 389 m 
14853. 800 m 
1 / 38, 185 



TABLE IX 

LATITUDE AND DEPARTURE CORRECTIONS 



Traverse lines 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 



X-correction (m) 

- 0. 031 

- 0. 035 

- 0.061 
- 0 . 010 

- 0. 069 

- 0. 090 

- 0. 022 
- 0. 022 

- 0. 029 

- 0. 032 

- 0. 018 

- 0. 069 

- 0. 084 

- 0. 090 



Y-correction (m) 

+ 0. 007 
+ 0. 008 
+ 0. 012 
+ 0. 002 
+ 0. 015 
+ 0. 019 
+ 0. 005 
+ 0. 005 
+ 0. 006 
+ 0. 007 
+ 0. 004 
+ 0. 015 
+ 0. 018 
+ 0 . 020 
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ADJUSTED COORDINATES BY COMPASS RULE 
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Coordinates of stations 1 and 15 were known and held fixed. 



The minimization of Equation 4.12 can be done by taking partial derivative with 
respected to each ilnknown variable (X) 

= d<t>/dx = 2 X T B T WB - 2F T WB = 0 (4.13) 

Transposing Equation 4.13 and rearranging yields 

(B T WB)X = B T WF (4.14) 



or 



NX + U = 0 (4.15) 

T T 

where N = (B WB) is the normal matrix of dimension u by u, and U is B WF. Then 

X = -N _1 U (4.16) 

For the adjustment in traverse, the vector V in equation 2.13 is represented the 
residual of observed angles and distances. If there are n observed angles in a traverse, 
there will be n - 1 observed distances and the number of residuals becomes 2n+ 1 
which includes the residual of angles and distances. 

There are two types of condition equations in the adjustment of observation 
equations: the angle condition and distance condition. 

From Figure 4.2 the angle condition can be expressed by 



v ia P; * (®ij * ®ik) (4-17) 

where vj a is a residual of observed angle, ay is a forward azimuth, a^, is a backward 
azimuth and Pj is an observed angle. Equation 4.17 is suitable when ay ^ a^. If ay 
< ajj., the quantities in parenthesis must be added by 360°. 

The Equation 4.17 can be expressed in coordinates of the two stations as 



„ 1 X. - X. , X. - X. 

v ia = Pi - ( tan * ^ TT ) - tan^M— J ) ] 

j ' 1 i 1 k * * i 



(4.18) 
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Figure 4.2 Determination of Angle and Distance from Coordinates. 



and 



Si = 



, X. - X. , X. - X. 

[ tan (— 1 — ) - tan (— L ) | 

Y. - V. Y, - Y. 



(4.19) 



J • k 

When linearized by using liquation 2.12, Liquation 4.19 yields 



gj = gjtXj', Yj', Xj', Yj', X k ', Y k ') + a,5Xj + a^Yj + a 3 5Xj + a 4 SYj + 
a 5 5X k + a 6 «Y k <«»> 

where Xj', Yj', Xj', Yj', X k ', and Y k ' arc the estimated station coordinates. By 
substituting the linearization of function gj, liquation 4.18 can be written as 



v ia + a l 6X i + a 2 5Y i + a 3 5X j + a 4 5Y j + a 5 6X k + a 6 5Y k = F 



(4.21) 



where 5Xj, 5Yj oX k , and SY k arc unknown parameters (correction in X and Y 

coordinates), cocfiiccints aj, a 2 a^, and a^ arc the partial derivative of function gj 

with respected to Xj. Yj X k , and Y^, respectively. 

For the distance condition, it can be expressed as 



V W - j; • I (Xj - Xj) 2 + (Yj.Yj) 2 ! 1 '' 2 



(4.22) 
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where 



hi = [ (Xj - Xi) 2 + (Yj - Yj ) 2 I 1 / 2 
The linerized form of Equation 4.23 is then given as 



(4.23) 





X k \ Y k ') + b,8Xi + b 2 8Yi + b 3 8X- + b 4 8Y- (4.24) 



Substituting the linearized form from Equation 4.24, Equation 4.22 becomes 



v id + bjSXj + b 2 8Yi + b 3 8X- + b 4 8Y- = F 2 (4.25) 

where FI and F2 represent the constant form for Equations 4.21 and 4.25. Table XI 
shows the partial derivatives for Equations 4.21 and 4.25. 

Sometimes, the station positions which is determined by intersection from 
traverse station is also to be adjusted simultaneously with traverse station positions. 

In this case only the angle conditions are added to observation equations of traverse. 

In the data adjusted under this thesis, there are 15 observed angles and 14 
oberved distances which consisted of 13 unknown points, consequently, there will be 29 
observation equations included 26 parameters of 8X and 5Y. Thus, they can be 
expressed in matrix form of Equation 2.13 as 



V 29,l + B 29,26 X 26,l “ F 29,l 



(4.26) 



where Vj j, v 2 j, v^ j are the residuals of angles; v^ j, Vj -j j, v^ j are the 

residuals of distances; B is the coefficeint matrix (consistings of a's and b's) of 
parameter X (Table XI); and F is the constant vector. 

When three intersection points with 10 observed angles were added (Figure 3.1) 
to adjust the coordinates, the observation equations have 39 equations including 32 
unknown parameters. 



V 39,l + B 39,32 X 32,l “ F 39,l 



(4.27) 
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TABLE XI 

THE COEFFICIENTS OF ANGLE AND DISTANCE CONDITIONS 



a 

a 

a 

a 

a 

a 

b 

b 

b 

b 



1 

2 

3 

4 

5 

6 
1 
2 

3 

4 



dX ± 

d *i 
5Y i 

dgj _ 

d<3i 

dY i = 
dgl 
5x k ” 

SY k 

dX ± 
dh^ 

dT ± 



dh.^ 



= + 



= + 



= + 



1 _ 



Y • - Y . 

13 11 + 

72 + 



( s i j ) 



= + 



= + 



< s ij) Z 
Y i - Y i 

( s ij) 2 
X i ~ X i 

Y k - Y i 
( s ik> z 
x k - x i 

( s ik) 2 
x i - x i 

s u 

Y i - Y i 

s ij 

x i - x i 

s ij 

Y i - Y i 
S ij 



Y k - Y k 

( s ik) 



X, - X, Xv, - Xi 



< s ik> 



By using Equation 4.16, N'* is the cofactor matrix and the diagonal terms of this 
matrix gives the variances of the adjusted coordinates. To obtain the residuals of all 
observed quantities, the reverse process must be done. The correction of X and Y 
adjusted coordinates and their standard deviation (<t's) are shown in Table XII. And 
the adjusted standard deviation of observed quantities were obtained by multiplying 
standard deviation of unit weight [(V ^ WV / r)*^] to a squares root of diagonal 
element of B(B T WB)* 1 B T matrix (Table XIII). 
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TABLE XII 

ADJUSTED COORDINATES BY OBSERVATION EQUATION 
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TABLE XIII 

ADJUSTED STANDARD DEVIATION OF ANGLES AND DISTANCES 



umber 


Angles 


Di stances 




seconds 


meter 




(± ) 


( 


1 


9. 8 


0. 029 


2 


8. 9 


0. 035 


3 


9. 1 


0. 045 


4 


9. 2 


0. 024 


5 


9. 3 


0. 048 


6 


9. 4 


0. 056 


7 


8. 4 


0. 030 


8 


8. 5 


0. 029 


9 


7. 8 


0. 032 


10 


7. 1 


0. 034 


11 


9. 6 


0. 028 


12 


9. 7 


0. 048 


13 


9. 8 


0. 054 


14 


6. 8 


0. 058 


15 


6. 8 
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V. ANALYSIS OF RESULTS 



A. COMPARISON OF ADJUSTED COORDINATES 

If the given coordinates of the control points are assumed to be error free, then 
the accuracy of the traverse station coordinates depends only on the accuracy of 
distance and angle measurements. The adjusted traverse coordinates obtained by the 
approximation method are of a lower order of accuracy as only the errors in the 
misclosure in azimuths and distances were determined. These errors were distributed 
by assuming that all observed quantities had an equal probable occurrence. 

The least squares adjustment method provides a better approximation of the true 
value. Therefore, the adjusted traverse coordinates obtained by this technique provided 
better estimates for position of all traverse stations and the accuracy of the adjustment 
can be checked and statistically tested. 

After the 13 adjusted traverse station (from stations 2 to 14) coordinates 
obtained by the approximation method and by the least squares method were 
compared, the difference in coordinates at each station were computed and plotted 
(Table XIV and Figure 5.1). The largest difference was at station 14. Because stations 
1 and 15 are held fixed, the least squares techniques adjusts simultaneously errors in 
azimuths and distances, while the adjustment by approximation adjusts errors in 
azimuths and distances sequentially. Consequently, the largest difference in traverse 
distances occurs at the last station (station 14) before closing of traverse at the fixed 
station. 

When the standard deviation of observed quantities before the adjustment 
(Tables III and IV) were compared with those obtained through adjustment (Table 
XIII), the standard deviation of all observed quantities in Table XIII showed 
increments. That means, the estimated standard deviations were optimistic. 

In this thesis, the three intersection points were also adjusted. The adjusted 
coordinates of these were compared to NOS results. The largest difference occurs at 
point no. 3 (5x = + 0.140 m; 5y = - 0.158 m). The standard deviation of adjusted 
coordinates at this point are crxj = ± 1.87, cry 3 = ± 0.96 m (Table XV). 
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TABLE XIV 



COMPARISON OF ADJUSTED COORDINATES/DISTANCES 
OBTAINED BY APPROXIMATION AND LEAST SQUARES METHODS 



Differences* 

Coordinates Distances 



stn. 




dx 




2 




f>%9 


_ 


3 


- 


0. 015 


- 


4 


+ 


0. 004 


- 


5 


- 


0. 013 


- 


6 


+ 


0. 002 


- 


7 


+ 


0. 008 


+ 


8 


+ 


0. 001 


- 


9 


+ 


0. 005 


- 


10 


- 


0. 008 


- 


11 


- 


0. 003 


- 


12 


- 


0. 014 


+ 


13 


- 


0. 005 


+ 


14 


- 


0. 019 


- 



* Approximation minus 



dy 

(m) 


( m) 


o; 007 


0. 020 


0. 017 


0. 023 


0. 022 


0. 022 


0. 012 


0. 018 


0. 008 


0. 008 


0. 024 


0. 025 


0. 005 


0. 005 


0. 013 


0. 014 


0. 023 


0. 024 


0. 034 


0. 034 


0. 003 


0. 014 


0. 041 


0. 041 


0. 101 


0. 103 


least 


squares solution 



B. ANALYSIS OF THE REFERENCE VARIANCE OF UNIT WEIGHT 

The weight matrix was set for the least squares adjustment by using Equation 
2.1 1. The <Tq (a priori reference variance of unit weight) was assumed as 1. The 
result of Gq = (V 1 WV / r) was obtained after adjustment. The value of Cq (a 
posteriori reference variance of unit weight) can be used to evaluate the weighting 
scheme used in the least squares adjustment. 

The standard deviation of the observed angles (<r a ) and distances (<r^) used in the 
least squares adjustments (solutions 1, 2, and 3) and the corresponding a posteriori 
obtained for these solutions are listed in Table XVI. As the a posteriori G () - for the 
solution no. 3 is closest to the assumed a priori <7 q (= 1), the weight (or the standard 
deviations) used for observed angles and distances in this case seem the most realistic. 
Further statistical testing for Gq done by Chi-squares or F-test was not earned out 
under this thesis. 
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DISTANCE DIFFERENCES (M) 

0.00 0.02 0.04 0.06 0.08 0.10 



Oi 

* 




Figure 5.1 Comparison of Adjusted Distances Obtained 
by Approximation and Least Squares Methods. 
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TABLE XV 

COMPARISON OF COORDINATES AT INTERSECTION POINTS 







X(m) 


crx( ±m) 


Y(m) 


cry( ±m) 


At 


Point No. 
NOS 

Adjusted 


1 

326616. 356 
326616. 470 


1. 72 


3071294. 683 
3071294. 635 


1. 03 




8x 


- 0. 114 


Sy 


+ 0. 048 




At 


Point No. 
NOS 

Adjusted 


2 

327056. 079 
327056. 183 


1. 68 


3070340. 474 
3070340. 444 


1. 14 




8x 


- 0. 104 


5y 


+ 0. 030 




At 


Point No. 
NOS 

Adjusted 


3 

325378. 101 
325377. 961 


1. 87 


3077263.378 
3077263. 536 


0. 96 




5x 


+ 0. 140 


5y 


- 0. 158 





TABLE XVI 

COM PARI SION OF VARIANCES OF UNIT WEIGHT 


Solutions 


Variance used in Adjustment 


V T WV/r 




2 


/N 2 




ff a a d ff 0 


*0 


1 


2" ( lOppm + 0. 5 cm) 1 


15. 08 


2 


5" (lOppm + 1. 0 cm) 1 


3. 98 


3 


10" ( lOppm + 2.0 cm) 1 


1. 32 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

By using the weight of observed quantities for the adjustment, the traverse 
station coordinates computed and adjusted by the least squares observation equations 
method were more accurate than those obtained by the approximation methods. Even 
though the observation equations method may require a greater number of equations 
than the condition equations method, the processing of the data for adjustment is 
easier and the corrections in X and Y coordinates can be directly obtained through 
iterative solution. This method is suitable when a computer with a memory capacity of 
over 500 K bytes is available. However, for local work or a relatively short traverse, an 
approximation method is commonly utilized when economic and logistic criteria are 
considered. 

B. RECOMMENDATION 

The INDTRA Fortran program written for this thesis is automated for handling 
only two kinds of survey techniques: traversing and intersection. With the computers 
available at NPS, the development of adjustment programs for covering a wide range 
of survey techniques should be done to use and continue analysis of the mixed kind of 
survey techniques including traverse, triangulation, trilateration, resection, and 
intersection. 
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APPENDIX A 
LINEARIZATIONS 

This is a linearized form of m functions in n unknowns. 



y l = r 1 (x 1 ,x 2 ,...,x n ) 

y 2 = f 2 (x 1 ,x 2 ,..,x n ) 



y m = U x l’ x 2’- x n) 



yO = 



y.° 




f^Xj 0 , x 2 °, ... 


x n °) 


>2° 




f 2 (x 2 °, X 2 °, .. 


- X n°) 


y m ° 




f (x X °, 

m v m ’ m ’ 


X n°> 




dY 

~dX 



8 _h 


d h 


5y, 


^Xj 


dx 2 


8x n 


dy 

J m 








dy 2 


<3x n 



AX = 



Ax. 

AxJ 



Ax, 



The general form of linearized functions becomes 
Y = Y° + J yx AX 
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APPENDIX B 

TRAVADJ FORTRAN PROGRAM 



This program is used for computing and adjusting the traverse station position by 
approximation method (Compass rule). 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



FORTRAN PROGRAM "TRAVADJ" 

THIS PROGRAM IS USED FOR 

1) REDUCE SLOPE DISTANCE TO ELLIPSOIDAL DISTANCE 

2) DETERMINE GRID DISTANCE 

3) COMPUTE CLOSE TRAVERSE 

4) ADJUST COORDINATE BY COMPASS RULE 



II a II 

2) VALUE OF 1/F (EX. 1/F = 294.978698 

3) CENTRAL SCALE FACTOR 

4) LAT° AT STARTING AND CLOSING POINT 

"idl.iml, sl.id2.1m2,s2" 

5) GRID N. AND £. OP 4 fcNOWN STATION 

"qnl ,gel ,qn2,qe2,qn3 ,ge3.gn4,gn4" 

6) NUMBER OF MEASURE DIStANCf 



7 ) 



) 



10 ) 

11 ) 

12 ) 



ELEVATION AT FIRST OCCUPIED POINT 

NAME S 0F ALL STATIONS 
INDICATOR VALUE 1 = NO VERTICAL ANGLE 
2 = VERTICAL ANGLE 

DIFFERENT IN ELEVATION BETWEEN TWO STATIONS 
SLOPE DISTANCE (WITH UNIT FEET OR METER) 

HORIZONTAL ANGLES 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

INPUT DATA 

1) SEMIMAJOR AXIS OF REFERENCE ELLIPSOID USED 



ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 
c 



ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 



c sfac = scale factor coorection 

c hgd = horizontal distance c 
c todis= total distance in traverse c 
c cegrid = grid east of traverse station c 
c cnqrid = grid north of traverse station c 
c difazi = azimuth misclosure c 
c difdis = distance misclosure c 
c coraz = angular correction per station c 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
C c 
c gridaz = subroutine for computing azimuth c 
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oooo nnnnnnnnnnnnnn 



utm 



dmsr 



rdms 



cccccccco 



between two traverse stations 
subroutine for computing the grid 
"Coordinates from known distance 
and azimuth 

subroutine for converting the angle 
from degrees, minutes, and seconds to 
radians 

subroutine for converting the angle 
from radians to degrees, minutes, and 
seconds 



ccccccccccccccccccccccccccccccccccccccco 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 

C 

C 



C 

c 



c 

c 



DIMENSION LAT(3) ,K(3) ,HGD(30) , HAGL( 30) , HANGf 30) 

DIMENSION DHAGL(3o) , mAAGL( 30 ) ’SHAgL( 30) ,NAMt(3,30) , NAME1 ( 5 ) 
DIMENSION NAME3( 2) ,f^lAME4( 5) ,MDlST( 30) 



* 

* 

* 

* 

* 

* 



VARIABLE DECLARATION 



DOUBLE PRECISION AN90R,AN180R,AN360R.AZC .AZFIR.AZIMUT .AZLAS , 
CEGRld( 30 T, CCiAZRC 30j .COftAZ ,CN6RID(30 j , 

distan,difAzi,deltax(30).d£ltay(30),dOmmyi, 

DUMMY2 DUMDIS DIFDIS.TODiS 



DOUBLE PRECISION 



DOUBLE PRECISION 



* 

* 

* 

* 



DOUBLE PRECISION 



NUM0 , NUM360; NUM1 . SUMDX J SUMDY^ STDD( 30 ) , STDA( 30) , 
XGD,VGD,NUMf80,NL)M90,Wf*RED(3;? 

ANGS( 30) ,CFAZS( 30) ,DUM1D .DUM1M. DUM1S 
,EAZIS,EXZ2S,C0RDxi,C0RD?l,TEM$,AZS21,AZS34 

SUMIX1 . SUMIX2 , SUMFY , SUNFX , SUMCOY . SUMCOX , 
SSTA,S$TD.SIG)i.SIGY,SIGXY’SEMAJ.$EMIN, 

SETA, SETAA.SSl6X.SSiGY.SSiGXY.DfGl, MINI. 

S ECOKl 1 , SEMAJ1 , SEMIN 1 , T!iMP01( 36) , TEMP02(30) , 
DDIFFX .DDIFFY 

REAL*8 A.C,D,H,R,V,X,Y,GD,DH,ZD.CON.CUV a UT,PHE,CUVRE.ELEV 

real*8 HSiST ’ s6l$T !m6i 6T , fSQft ,GftAMMA,AN6 ,AN(jV , Pi ,AGU , ftl , F2 

REAL*8 K,KO.fl.E2.0l,Q2,ftR,VV ECEN.FAC.SfAC.H6D 

REAL*8 F CH6D.fiSDiST DDM.HA.HAGL.CADISt.DG.fiMS.DS.DDANG.MMANG 

REAL*8 GNl.GNf.GN3,GN4,Gfl,6E2.Gf3.GE4.AZ2i,AZ34.6lST21,DIST34 

REAL*8 Sl.62.Dl.D2 SV.^ANG SHA6L,SANG,6D,MM,SS,HANG,RRR 

REAL*8 DUM1,&UM2,DUM3 DUM4 XXXX.YYYY 

REAL*8 IDD1 IMM1 , IDD2 IMM2,DVV,MVV, CUVRER ,CUVR, CURVR 

REAL*8 ZINEjKOSE 

INTEGER I,N.ID1.ID2.IM1.IM2.DV,MV.DHAGL.MHAGL,DANG.MANG,UNIT 
INTEGER ANG6(30i.AN6M(3fl).A6Jl’ERft,AZD2i.AZM2i.AZD34,AZM34 
INTEGER SSHAGLX36).CFAZDjr60LcfAZM(30J.MiN2,DE62, 

EAZ1D, EAZ lM , EAZ2 D , EA^2M , TEMD , Tf MM 

DATA DUM1/0. 30480061D0/ .DUM2/500000. 0D0/,DUM3/0. 000001D0/ 

DATA DUM4/3600. 0D0/ , SIGANG/5. 0D0/.DUM5/0. 000005D0/,DUM6/0. 005D0/ 
DATA TODIS/O. 0D0/,DUM10/0. 0000048481368D0/ 



DATA NUMO/O. 0D0/.NUM360/360. 0D0/.NUM1/1. 0D0/.SUMDX/0. 0D0/, 
* SUMDY/0. ODO/ , NUM180/180. 0D0/.NUM90/90. OdO/ 



READ( 5,10 
READC5 15 
READC5 15 
READC5 16 
READ(5,16 
READC5 16 
READ(5 16 



A a F,K 0 
l6l IM1.S1 
ID2, IM2 ,S2 
GN1.GE1 
GN2 GE2 
GN3 ,GE3 
GN4 GE4 
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READ( 5 , 1 1 )N , ELEV 

CALL GRIDAZ (GE1 , GN1 , GE2 , GN2 , AZ21 , DIST21 ) 
CALL GRIDAZ (GE3 |GN3 ,GE4 ,GN4 ,AZ34 ,DIST34) 

IDD1 = FLOAT C ID1 ) 

IMM1 = FLOAT (IM1) 

CALL DMSR ( IDD1 . IMM1 ,S1 ,D1) 

LATfl) = 6 1 
IDD2 = FLOAT! I D2) 

I MM2 = FLOAT U M2) 

CALL DMSR !IDD2.IMM2,S2,D2) 

LAT ( 2 ) =62 

LAT { 3 ) = (Dl+D2)/2. 0D0 
ESQR = 2.0D0*(1.0D07 fJ-( 1.0D0/F)**2 
X = A*DSQRT!l. ODO-ESQR) 

Y = 1.0D0-(ESQR*(DSIN(LAT(3)))**2) 

R - X/Y 



DETERMINATION OF THE SCALE FACTOR FOR UTM. 



El = DABS(DUM2-GE2' 
E2 = DABS(DUM2-GE3; 



ECEN = ESQR/C1. ODO-ESQR) 



)1 = DUM3*E1 
2 = DUM3*E2 

PRIME = ((Ql**2)+(Ql*Q2)+(Q2**2))/3. 0D0 



DOl M=1 ,3 

RR=A/( 1. 0D0-ESQR*(DSIN( LAT(M)))**2)**0. 5 

FI = ( 1. ODO+ECEN*DCOS! LAT(M) ) )*( 10. 0D0**12) 
F2 = 2. ODO*(RR**2)*(KO**2) 

FAC= F1/F2 



CON T K ffif K0 * (1 ' 0D0+FAC *Q pRIME+ (°- 00003D0*(QPRIME**2))) 



SFAC = K(3) 

SFAC = 6. UD0/((1. 0D0/K(l))+(4. 0D0/K(3) )+( 1. 0D0/K(2))) 



READ(5 ,20) UNIT 
READ(5!l8)(NAMEl(L),L=l,5) 

DETERMINATION OF THE HORIZONTAL DISTANCES 

DO 1000 J=1,N 

READ! 5 , 14)!NAME! L, J) , L=1 ,3) 

READ! 5 , 12 )I ,DH,DV MV SV 
READ!5.13) $DI$T 
MDIST(J) = SDIST 

IN CASE OF THE LENGTH'S UNIT IS IN FEET, THEN CONVERSE TO METER 



I F( UNIT. EQ. 1) THEN 
SDIST = SDIST*DUM1 
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DH = DH*DUM1 
ELEV = ELEV*DUM1 
END 4F 



READ( 5^1 7 )DHAGL^J LMHAG^J ) , SHAGL(J) 
MM = DFLOATCMHAGlIj)) 



SS = SHAGL(J) 

CALL DMSR ( DD,MM,SS,RRR ) 
HANGfJ) = RAR 
STDA(J) = SIGANG 



DETERMINATION OF HORIZONTAL DISTANCES WHEN THE DIFFERENT 
ELEVATION IS APPROXIMATELY KNOWN. 

I F( I. EQ. 1)THEN 

IFfDH.NE. O.OJTHEN 
DDM = SDIST 
IFQDDM.GE. 3300.0)THEN 
DO 100 KK=1 ,3 



IN 



CURV =[0. 016192D0*DDM) 

CALL DMSR (NUMO,NUMO A CURV,CURVR) 
D=DSQRT(DDM**2-{DH*DCOS(CURVR)r*2)-DH*DSIN(CURVR) 



100 



CALL DMSR (NUMO.NUMO 
D=DSQRT 
DDM = D 

CONTINUE 

DDM 

-DH**2) 



HDIST = 

ELSE 

HDIST = DSQRT(DDM**2- 
END IF 



ELSE 

HDIST = SDIST 
END IF 
END IF 

DETERMINATION OF HORIZONTAL DISTANCE WHEN ZENITH DISTANCE IS KNOWN 

IF(I. EQ. 2)THEN 

DVV = FLOAT(DV) 

MW = FLOATfMV) 

CALL DMSR (DVWMVV,SV,VANG) 

CALL DMSR (NUM§0 .NUMO . NUM0 a AN90R) 

CALL DMSR ( NUM 18(5 , NUM6 , NUM6 ,AN180R) 

ZD = AN90R-VANG 
DH = SDIST*DSIN(ZD) 

DS = SDIST*DCOS[ZD) 

I F (S D I S T . GE. 3300. 01 THEN 
CUVRE = 0. 01392500*SDIST 
CUV = 0. 016192D0*SDIST 
CALL DMSR (NUMO , NUMO , CUVRE , CUVRER) 

CALL DMSR (NUMO’NUMO,CUV,CUVR) 

C 

GRAMMA 
HDIST 
ELSE 

HDIST = DS 
END IF 
END IF 



= (AN96r-ZD)+CU^RER 



(AN90R+CUVR) 

AN180R-(C+PHE] / s 

( SD I ST*DSIN( GRAMMA) )/DSIN(C) 
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HH= (ELEV+DH) 

determination^ the horizontal distance on geoid 

GD = R*HDIST/( R+H) 

DETERMINATION OF THE GRID DISTANCES 
HGD(J) = GD*SFAC 



TODIS = TODIS+HGDfJ) 

STDD(J) = HGD(J)*DUM5+DUM6 

WRITE(6,*) 'GRID DISTANCE =' ,HGD(J) 

WRITE( 6 ’ * ; 1 Y 

ELEV = HH 

1000 CONTINUE 

READ( 5 , 19 ) ( NAME3( L) ,L=1 ,2) , DANG , MANG , SANG 
READ( 5 ’ 18)(NAME4( L) ) L=1 | 5 j 

CALL OUTPUT! N, AZ21 ,AZ34 .GN2 ,GE2 ,GN3 .GE3 ,HGD,DHAGL.MHAGL,SHAGL, 
*DANG, MANG , SANd , NAMEl . NAME3 , AME4 , NAME ,MdlST , GN1 , GEl , GN4 , GE4 ) 
dDANG = DFEOAT(D4NG) 

MMANG = D FLOAT! MANG) 

CALL DMSR (DDANG, MMANG, SANG, RRR ) 

NN = N+l 
HANG(NN) = RRR 
STDA(NN) = SIGANG 



C 



C 



C 

C 

C 

C 

C 

c 



CALL 

CALL 



DMSR 

DMSR 



NUM360 ,NUM0 ,NUM0 ,AN360R ) 
NUM180 , NUMO ,NUM0 ,AN180R ) 



DO 2000 ADJ1 = 1,2 
XGD = GE2 
YGD = GN2 
AZFIR = AZ21 
DO 200 I = 1 ,N 

AZIMUT = AZFIR+HANG( I) 

IF (AZIMUT. GE.AN360R) THEN 
AZIMUT = AZIMUT-AN360R 
END IF 

IF(AZIMUT. GE. AN180R) THEN 
AZIMUT = AZIMUT-AN180R 
ELSE 

AZIMUT = AZIMUT+AN180R 
END IF 

COAZR(I) = AZIMUT 



CALL RDMS f AZIMUT,DUM1D,DUM1M,DUM1S ) 
CFAZD(I) = idlNT! DUM1D) 

CFAZM!li = IDINT(DUMIM) 

CFAZS(I) = DUM1S 



DISTAN = HGD! I) 

CALL UTM ( XGD, YGD, DISTAN, AZIMUT, DUMMY1 ,DUMMY2 ) 
WRITE(6,*)'I =' ,1 
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WRITEC6 ,*) ' GRID E=' . DUMMY 1. 
DELTAX(I) = DUMMYl-XGd 



GRID N=‘ ,DUMMY2 



_ DELTAYCI) = DUMMY2-YGD 

17RITE(6 *) DEPARTure =' DELTAX(I) , 1 LATITUde =' 
ZlNE = DSIN(AZIMUt) 

KOSE = DCOS(AZIMUT) 



SUMDX = SUMDX+DELTAX( I ) 

SUMDY = SUMDY+DELTAY( I ) 

CEGRID(I) = DUMMY1 
CNGRID(I) = DUMMY2 

CALL GRIDAZ ( XGD , YGD , DUMMY 1 , DUMMY2 , AZC , DUMDIS ) 

CALL RDMS ( AZC.DUMID^UMIM^DUMIS ) 

CBAZD(I) = IDINT(DUMID) 

CBAZM(I) = IDINT(DUMIM) 

CBAZS(I) = DUM1S 



IF ( ADJ . EQ. 1) THEN 



END I 



EMPOl(I) = CEGRID(I) 

TEMP02(I) = CNGRID(I) 



AZFIR = AZC 
XGD = DUMMY 1 
YGD = DUMMY2 
200 CONTINUE 

IF (ADJ1.NE. 1) THEN 

CORDX1 = ( SUMDX' 
CORDY1 = l SUMDY 
XGD = GE2 
YGD = GN2 



-(GE3-GE2) J/TODIS 
-(GN3-GN2) )/T0DIS 



DO 230 1=1, N 

DISTAN = HGD(I) 

XXXX = C0RDX1*HGD( I) 

YYYY = CORDYl*HGDa) 

DELTAX(I) = DELTAX(I) 
DELTAY(I) = DELTAY(I) 



XXXX 

YYYY 



CEGRID(I) 

CNGRID(I) 

DDIFFX 

DDIFFY 



XGD+ DELTAXfl' 

YGD+ DELTAYf I ' 
CEGRID(I)-TEMP01(I) 
CNGRIDC I )-TEMP02( I ) 



C 

c 



230 



XGD 

YGD 

CONTINUE 



= CEGRIDf I) 
= CNGRID( I) 



END IF 

AZLAS = AZFIR+HANG(NN) 

IF ( AZLAS.GE. AN360R ) THEN 
AZLAS = AZLAS-AN360R 
END IF 

IF (AZLAS. GT.AN180R) THEN 



, DELTAY ( I ) 
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c 

c 



AZLAS-AN180R 

AZLAS+AN180R 



AZLAS 
ELSE 
AZLAS 
END IF 

CALL RDMS ( AZLAS ,DUM1D.DUM1M,DUM1S ) 



TEMD = IDINt(DUMlD) 
TEMM = IDINT(DUMIM) 
TEMS = DUM1S 



DIFAZI = AZLAS-AZ34 
WPRED(l) = DIFAZI , N 

CORAZ = DI FAZ I/DFLOAT ( NN ) 

CALL RDMS ( DIFAZI ,DUM1D,DUM1M,DUM1S ) 
EAZ1D = IDINt(DABS( DUMlDn 
EAZ1M = IDINTC DABS(DUMIM)) 

EAZ1S = DABS(DUMIS) 

CALL RDMS ( CORAZ .DUM1D.DUM1M.DUM1S ) 
EAZ2D = IDINT(DAB$(DUM1D)) 

EAZ2M = IDINTCDABS(DUMIM)) 

EAZ2S = DABS(DUMIS) 



IFfADJl. EC 
WPREDC ‘ 
WPRED 
DIFDI. 
ERR 
END IF 



1)THEN 

= SUMDX - ( GE3-GE2 ) 

= SUMDY - I GN3-GN2 ) 

= DSQRTf WPRED(2)**2 + WPRED(3)**2 
= IDINT( TODIS / DIFDIS ) 



PRINT RESULTS OF TRAVERSE COMPUTATION 



CALL RDMS (AZ21 ,DUM1D.DUM1M,DUM1S) 
AZD21 = iDINTtfiUMID) 

AZM21 = IDINT(DUMIM) 

AZS21 = DUM1S 

CALL RDMS (AZ34,DUM1D a DUM 1M,DUM1S) 
AZD34 = iDINT(d)UMlD) 

AZM34 = IDINT(DUMIM) 

AZS34 = DUM1S 



IF (ADJ1. EQ. 1) 
WRITEC8 ,29) 
WRITE(8 27^ 
WRITE* 8*30 
WRITE(8 31 
ELSE 

WRITE(8,28 
WRITE(8 27 
WRITE* 8,30 
WRITE(8 40 
END IF 
WRITEC8 , 32 
WRITE* 8 , 33 
WRITE(8 30 
WRITE* 8,34 
WRITE* 8 ’ 35 



THEN 



AZD21 ,AZM21 ,AZS21 
GE2 , G^I2 



DO 240 1=1, N 
II = 1+1 

CALL RDMS(HANG( I) ,DUM1D,DUM1M.DUM1S) 
ANGD(I) = IDINT(DUMID) 

ANGM(I) = IDINT(DUMIM) 

ANGS(I) = DUM1S 

SSHAGL(i) = IDI NT ( SHAGL( I ) ) 
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IF (ADJ1. EQ. 1) THEN 
wRITE(8,36‘ 

ELSE 



cTO!iTOM!f L(I) ' CFAZD(I) ’ 



C 

c 



240 



WRITE(8,37^ANGpn)^ANGMfI) ,ANGS(I) , CFAZD( I ) , CFAZM( I ) , 

END IF ( { ; 

WRITE(8 , 38) II ,CEGRID( I) ,CNGRID( I) 

CONTINUE 



CALL RDMS ( HANG(NN),DUM1D,DUM1M,DUM1S) 
IF (ADJ1.EQ. 1) THEKl 



DHAGL( 

MHAGLf 



= DANG 
= MANG 



ELSE 



SSHAGL(NN) = IDINT(SANG) 

_WRITE(8,36)DHAGL(NN),MHAGL(NN),SSHAGL(NN) JEMD JEMM JEMS 
CALL RDMS (HANG(NN) ,DUM1D,DUM1M,DUM1S) 



END IF 



ANGD(NN) = IOINTIDUMID, 

ANGM(NN) = IDINT(DUMIM) 

ANGSfNNi = DUM1S 
WRITE(8,37)ANGD(NN),ANGM(NN),ANGS(NN) JEMD JEMM JEMS 



3.1) THEN 
i , 4 1 ) EAZ ID , EAZ 1M , EAZ 1 S 



C 

C, 

c‘ 

c 



WRITE(8,30) 

WRITE(8,39;AZp34,AZM34,AZS34,GE3,GN3 
WRITEC8 *Y ' y 
IF (ADJl. EQ. i: 

WRITE(8, , . , 

WR ITE( 8 , 42 ) EAZ 2D , EAZ2M , EAZ2S 
WRITE(8 43JTODIS 
WRITEC8 44)DIFDIS 
WRITE(8,45) ERR 
END IF 

CORRECTED OBSERVED ANGLES 

DO 250 1=1, NN 

HANG(I) = HANG(I)-CORAZ 
250 CONTINUE 

SUMDX = 0. 0D0 
SUMDY = 0. 0D0 



2000 CONTINUE 



SUMDX = GE3-CEGRID(N 
SUMDY = GN3-CNGRIDCN 
DIFDIS= DSQRT ( SUMD 



WRITEf 8,41) EAZ1D. 
WRITE(8’44; DIFDlS 



**2 + SUMDY**2 ) 
EAZ1M, EAZ1S 



10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 



FORMATS 5X.F15. 5 . 5X . F15, 10 , 15X , F12. 8) 
FORMATS 26* ,I4,3$X.f7. 3) ' 

FORMATS 20X Il’4X.f9. 3 , 11X ,213, F10. 4) 
FORMATf 25X! Fl$. 6) 

FORMATS 6X.3A4) 

F0RMAT(5X!2I3J14.8) 

FORMATmO#,F1!d. 5.15X.F15. 5) 

FORMAT ( 25X*2I3 , FiO. 4) 

FORMAT f 5A4) 

FORMAT(2A4,17X,2I3,F10. 4) 

FORMAT T 1 1 ) 



' 27 FORMAT( 16X , 1 ::::::::::::::::::::::::::::: 1 ,///) 
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28 

29 

30 

31 
40 

32 

33 

34 

35 

36 

37 

38 
C39 
C41 
C 

C42 

C 

C43 

C44 

C45 



FORMAT(16X, 'ADJUSTED TRAVERSE COMPUTATION 1 ) 
FORMAT? 15X 'UNADJUSTED TRAVERSE COMPUTATION') 

cnDMATr i 



FORMAT? 

* I N 

FORMAT? 1 j 5TN 
MAT? - ' j ^TN 



rur\r 
*. ) 
FOfo 
*TEi 
FOR 



flAT 



( 



FORMAT? 1 ; 
*N0RTH : 1 
'I 



0. 



OBSERVED 
CORR. OBS. 
HOR. ANG. 

D M S 



: FORWARD 


: GRID : 


: FORWARD 


: GRID : 


: AZIMUTHS 


: DIST. : 


: D M S 


: (M. ) : 



COORDINATE (M 
ADJUSTED COORDINA 

GRID EAST : GRID 



I 1 

I I ,3X 



) 13X,I3 A 1X,I2,1X.F5.2.34X. . 
' 36X! Fid. 3!lX Fll.3.1*.'l') 

, r : Y ,i3,lX,l2,iX,I2,4x,i3,lX, 



'I') 

12 , IX , F5. 2,1X,F8. 3,25X 



FORMAT 
FORMAT, 

FORMAT? 1 
* 1 1' ) 

*25X A ( ' 1 ' ’ 3X > ' : ' ,13, IX, 12, IX, F5. 2, IX, 13, IX, 12, IX, F5. 2.1X.F8. 3, 
FORMAT? 1 1 ' 1 3 ' : ' 36X, F10. 3 . IX , FI 1. 3 , IX . ' I ' ) 

FORMAThx/KNdWN dATA' , 7X , li , li , 12 . 1>? . F$. 2 , 11X , F10. 3 , IX , Fll. 3.//) 
^FORMAJfj ANGULAR ERROR =' , 2x , 13 , IX , ' D' | IX , 12 | IX J ' M 1 , lft , F5. 

f6rmAt( ' At 

*2 a 1X. S 

f6rmAt? 

FORMAT? 

FORMAT? 

STOP 
END 



l 



NGULAR CORR./ STATION =' ,2X , 13 , IX , ' D',1X , 12 , IX , ' M',1X , F5. 



OTAL DISTANCES 
LINEAR ERROR 
LINEAR CLOSURE 



=' ,2X,F10. 3) 

=' 2X F10. 3) 

=' 3X ' 1 r 16,///////////) 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

C SUBROUTINE TO CALCULATE GRID AZIMUTH AND DISTANCE BETWEEN C 
C TWO KNOWN STATIONS WITH THEIRS GRID NORTH AND GRID EAST. C 

c c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

SUBROUTINE GRIDAZ (XI, Y1,X2,Y2, ANGLR, DIST12) 

c 

c 

REAL*8 XI ,X2 . Y1 , Y2 .ANGLR , ANGO ,ANG90 ,ANG180 .ANG270 , DIFX , DIFY .DIST12 
REAL*8 ANd90ft,AN18dR ,AN2^0R 



C 

c 

c 

c 



DATA ANGO/O. ODO/ .ANG90/90. 0D0/.ANG180/180. ODO/ .ANG270/270. ODO/ 



CALL DMSR 
CALL DMSR 
CALL DMSR 

DIFX = X2-X1 
DIFY = Y2-Y1 



ANG90 , ANGO .ANGO .ANG90R ) 
ANG18d ,ANGd ,ANGd ,AN180R 
ANG270, ANGO ANGO, AN270R 



C 

C 



IF? (DIFX. EQ. ANGO). AND. (DIFY. EQ. ANGO)) THEN 
ANGLR = ANGO 

ELSE IFfDIFX. EQ.ANGO) THEN 
IF(DIFY.GT.ANGO) THEN 
ANGLR = ANGO 
END IF 

ANGLR = AN180R 

ELSE 

IF(DIFX.GT. ANGO) THEN 

ELS^GLR = ANG90R-DATAN(DIFY/DIFX) 

ANGLR = AN270R-DATAN(DIFY/DIFX) 
END IF 

END IF 

DIST12 = DSQRT(DIFX**2+DIFY**2) 
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RETURN 

END 



0000000000000000000000000000000000000000000000000000000000 
0 (3 

0 SUBROUTINE OUTPUT TO PRINT OUT THE RESULT FOR USING IN 0 

0 THE NEXT COMPUTATION. 0 

0 0 

(3 (3 (3 (3(3 (3(3 (3 (3 (3(3 (3 (3 (3(3(3 (3 (3(3 (3(30(3 (3(90(3(3 (3(3 (3(3(3(3 (30(3(3(3 (3@P(3(3(3@P(3(9(3@(3(3(3(3@P@ 



c 

c 



c 

c 



10 



c 

c 



SUBROUTINE OUTPUT ( N . A , B , G2 , E2 . G3 . E3 , RD , DH , MH , SH , DA , MA , SA , NAME1 , 

*NAME3,NAME4,NAME,MDfSf,(al,£l,<a4,£4) 



im 

G,G2.G 



DH(N),MH(N),SH(N) ,NAME(3,N) , NAME 1 ( 5 ) , NAME3 ( 2 ) 



DIMENSION RD 
\NAME4(5).MD: . „ 

kEAL*8 V A\.&,DG,G2 a G 3.E2.E3 .RD.AA.BB.HH.SH ,SA,SSA,SIGA, SIGD 
REAL*8 Ml!) I $T , MD I § , R&D , (-il , 64 , £ 1 , £4 , ClUMl 
INTEGER N,DH,MH,DA,MA 

DATA SIGA/5. 0D0/.DUM1/0. 30480061D0/ 



WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 



= fl+i 



,240; 

250 

260 

250 

270 

280 

290 

280 

240 

,280 

300 

,310 

310 

300 

280 

240 



'NAME 

'NAME 

'NAME 

NAME4 






WR 

WR 

WR 

WR 

WR 

WR 

WR 

WR 

WR 

WR 

WR 



TE( 

TE( 

TE( 

TE< 

TE 

TE 

TE 

TE? 

TE( 

TE 

TE 



,230' 
,100 
,110 
120 
110 
125 
,140 
130 
140 
,100 
,140 



N1 


= N+l 


DO 


10 


C\J 

II 

l — l 




WRITE! 6 , 


* 

CONI 




r 1 -- 


WR] 


[TE( 


'6,150 


WR] 


[TE( 


6 150 


WR] 


[TE( 


6 140 


WR] 


[TE( 


6 100 



1) 
)(NAME( 



^t5!^:^!(^4 E a^= L i:d!bW!sr (N) ' sH(N) 



WRITE(6,230) 
WRITE(9 1 170) 



61 



c 

20 

C 

C 

40 

100 

110 

120 

125 

130 

140 

150 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 



WRT 

WRF 

WRI" 

WRF 

WRF 

WRF 

WRF 

WRF 

WRF 



'9,180' 
'9,490 
'9,180 
'9 170 
'9 200 
'9 210 
'9 200 
'9 170 



E(9,200; 



DO 20 J=1,N-1 

MDIS = MDIST! J)*DUM1 

RRD = RD! J)/DUM1 

SIGD = (RD(J)*0. 00001D0)+0. 01D0 

* ; 2 R 2 D? iSOT . 2 ’ , ( N A ME ( L , ) , L= ! , 2 ) , MDI S , 

CONTINUE 

MDIS = MDIST(N)*DUM1 

RRD = RD!N)/DUM1 x % 

^WRITE(9^220)(NAME(L,N),L=1,2),(NAME3(L),L=1,2),MDIS,MDIST(N) 

$IGD = (RD!N)*0. 00001D0)+0. 01D0 
WRITE(9,200) 

WRITE(9 , 170) 



FORMAT(2I3 . F6. 2) 

FORMATi 9X f********* 

**** *** T \ ’ 

FORMAT (^X, 1 * 

FORMAT^,'* MEASURED HORIZONTAL ANGLE * 



********************************************** 

* * * 

* * 



FORMAT! 9)! ******************************** DEGREES * MINUTES * SEC 
*ONDS * f r 

FORMAT^X,'* AT * TO 

FORMAT^!,'* * * * * 

FORMAT ! 9^ . 1 * 1 , 3X . 2A4 , 3X , 1 * 1 ,3X,2A4,3X, ,3X,I3,3X, ,4X,I2,3X, 
*'*' .2X.FS. 2,CX, f *') 

formAtc Sx . 1 kttkmhmmmmMmttmmmmmmmmm 

FORMAT! 5X,‘# » u 

* , . #') 

FORMAT! 5X'# DISTANCES 

*D DISTANCES # r ) 

FORMAT! 5X/# , # # 

* # # ) 

FORMAT! 5X '# FROM # TO # 

* .# - Wl... J\) 



# 

# 



MEASURED DISTANCES 
# 

METERS # FEET 



# 

» 

# 



REDUCE 

METERS 



FORMAT! 5X .2A4: 2A4, '# ' ,2X , F8. 3 , 2X , '# 1 ,2X , F8. 2 , 2X , 1 § ' , 2X , 
*F8. 3 , IX , # , £x , F&. 2 ,CX , A ) 

FORMAt('i') 

FORMAT( 9X I 

FORMAT! 9X, '* 

FORMAT! 9X, 1 * NAME 

FORMAT! 9X,'* OF 

FORMAT! 9X, '* 

FORMAT! 9X, '* KNOWN POSITIONS * GRID NORTH 
FORMAT( 9X , ' * 1 ,5A4 . 2X . FI 1. 3 ,2X, ' * ' 2X . Fll. 3 , 2X, ' * 1 ) 

RETURN 



* * l 

* COORDINATES *' 

********************************* I 

* * * I 

* GRID EAST *' 
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END 



(3(9(9(3(3(3(3(3(3(3Q(3(3PP(3(3(3(3(3(3(3(3(3(B(9(3(3PP(3(BPP 
0 0 

0 CONVERT DEGREES, MINUTES, 0 

0 AND SECONDS TO &ADIANS 0 

0 0 

0000000000000000000000000000000000 



SUBROUTINE DMSR ( D1,M1,S1,R1 ) 

REAL*8 D1,M1,S1,R1,PI1 ,N60 ,N3600 ,N180 

DATA N60/60. 0D0/.N3600/3600. 0D0/,N180/180. 0D0/ 

PI1 = 4. 0D0*DATAN( 1. 0D0) 

R1 = ((D1+(M1/N60)+(S17N3600))*PI1)/N180 

RETURN 

END 



C 

C 

c 



c 

c 



000000000000000000000000000000000 
0 0 

0 COMPUTE THE COORDINATES 0 

0 FROM ONE KNOWN POINT WITH 0 

0 MEASURED DISTANCE AND ANGLE 0 

0 0 

000000000000000000000000000000000 

SUBROUTINE UTM (XI , Y1 , DIS12 .ANG12R, X2 , Y2 ) 

REAL*8 XI ,Y1 ,X2,Y2 ,DIS12,ANG12R,DIFX ,DIFY 

DIFX = DIS12*(DSIN(ANG12R) ) 

DIFY = DIS12*(DCOS(ANG12R)) 

X2 = Xl+DIFX 
Y2 = Yl+DIFY 

RETURN 

END 

000000000000000000000000000000000 
0 0 

0 CONVERT ANGLE FROM RADIAN 0 

0 TO DEGREE, MINUTES, SECONDS 0 

0 0 

000000000000000000000000000000000 

SUBROUTINE RDMS ( RR1,D2,M2,S2 ) 

REAL*8 RR1 ,D2 ,M2 ,S2 , PI1 , N60 , N180 ,TDEG,DIF,TMIN 

DATA N60/60. 0D0/,N180/180. 0D0/ 

PI1 = 4. 0D0*DATAN(1. 0D0) 

TDEG= f N180*RR1)/PI1 
D2 = DFLOAT( IDINT(TDEG)) 

DIF = TDEG-D2 
TMIN= DI F*N60 
M2 = DFLOATflDINT(TMIN)) 

DIF = TMIN-M2 
S2 = DIF*N60 

RETURN 

END 
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APPENDIX C 

INDTRA FORTRAN PROGRAM 



INDTRA Fortran program is used for computing and adjusting traverse station 
position by least squares method of observation equations (indirect observation) and 
simultaneously adjust the intersection points which are observed at each traverse 
station. 



(3(3 (3(3 (3 (3 (3 (3 (3 (3(3 (3 @ (3 (3 @ (B (3(3@(3 (3(3@(3 (3 (3 (3 (30(3 
00 (30 

00 LEAST SQUARE ADJUSTMENT OF 00 

00 OBSERVATION EQUATIONS METHOD 00 

00 00 

000000000000000000000000000000000000000 



000000000000000000000000000000000000000000000000000000 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

N ,NA = NUMBER OF OBS. ANGLES 0 

Nftl = NUMBER OF RESECTION PT. 0 

ND = OBS. DISTANCES 0 

NT = OBS. DIST + OBS. ANG. 0 

NC = TRAV. STN. + RESECT. PT. 0 

(NOT INCLUDE KNOWN STN. ) 0 

NZ1 = NUMBER OF RESECT. ANGLE AT EACH TRAV. STN. 0 
EGRID = GRID EAST OF KNOWN POINT 0 

NGRID = GRID NORTH OF KNOWN POINT 0 

D.M a S = READ ANGLE IN DEGREES, MINUTES, SECONDS 0 
RAN6. DIST = GRID DISTANCES 0 

NCO, £C0 = GRID NORTH AND EAST OF INTERSECTION PT. 0 
STDA = STANDARD DEVIATION OF ANG. 0 

STDD = STANDARD DEVIATION OF DIST. 0 

NP2 = NUMBER OF UNKNOWN DX , DY 0 

CEGRID.CNGRID = TRAV. STN. CfiORDI. 0 

INDLSQ = SUBROUTINE ADJUSTING THE COORDINATES 0 

BY OBSERVATION EQUATIONS 0 

CALFM = SUBROUTINE FOR COMPUTING CONSTANT VECTOR 0 
CALAM = SUBROUTINE FOR COMPUTING THE COEFFICIENTS 0 



OF UNKNOWN PARAMETERS 



0 

0 



( 3(3 (3 (3 (3 ( 3(3 ( 3(3 ( 3(3 (3 (3 (3 (3 ( 3(3 ( 3(3 (3 (3 (3 (3 (3 (3 (3 (3 (3 ( 3(3 ( 3(3 (3 (3 (3 (3 (3 (3 (3 (3 (3 (3 (3 ( 3(3 ( 3(3 ( 3(3 (3 (3 ( 3 ( 3(3 



* 

* 

* 

* 

* 

* 

I 

* 

* 

* 

[ 

* 

* 



DOUBLE PRECISION 



DOUBLE PRECISION 



DOUBLE PRECISION 



NGRID(4).EGRID(4J.AZ21 .DIST21.AZ34, 
DIST34,XdD,YGD>Z[ f IR A NL)M360,NL)M0. 

AN360R ! DD ,MM , S$ , R , AN6l(40) . $TDA(40) . 

DUM10 , &AN6,DiST(4d) ,NCu,Ec6.CNGRID(40) , 
CEGRlflUO) ,STDD(40) ,AZIMUT,DlS(40) . 
coazr(4o) j D ummy 1.DUMMY2 ,azC,cfazs(4o) , 
DUM1D ,DUMiM,DUMl$ v DELX(40) ,DELY(4u 1 
SUMDX , SUMDY, CBAZSj 40 ) , TTCOk ( 40 j , TTCOYjf 40 ) , 
AZLAS.DIFAZf ,WPRED(3) . EAZISXORAZ . DIFDIS , 
EAZ2S ,T0DIS ,ASNY(40) ,ASEX(4D) , DUMi 1 , 

DUM12 

WMf 80 1 , 80 ). STD0(80) , FMf 80 ) . BMC 80 ' 80) . 

BMT(8d ,80 j ,BMTWM(od.8U) ,NM(80 A 8d) ,TM(80) , 
DELTA(DO) ,NMI(80,80) ,WKlO(428D) ,VM( 1 ,80) , 
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c 

c 



STDAD,XCOX(80) ,YCOY(80) 



c 

c 



c 

c 

ce@ 

coo 



cooo 

c 

c 

coooo 

c 

c 

c 



c 

c 

c 

c 



c 

c 



c 

c 



c 



* 



INTEGER 

* 

ye 

* 




DATA NUM360/360. 0D0/,NUM0/0. 0D0/ .DUM10/2. 0D0/ 
DATA DUM11/0. 00001D0ADUM12/0. 0260/ , SUMDX/0. 0D0/ 
DATA SUMDY/0. ODO/,TODiS/O.ODO/ 

DATA PW/0/ , PRINTD/1/ 

PW = 1 UNWEIGHT 

PRINTD = 0 NOT PRINT DETAILS OF MATRIX 



DETERMINE BACK AZIMUT AT STARTING STATION AND FORWARD 
AZIMUTH AT CLOSING STATION 



DO 1 1=1,4 




READ! 5,*) 


i NGRIDl 


READ(5 * 


i EGRIDI 


CONTINUE 




CALL GRIDAZ i 


f EGRIDC 


CALL GRIDAZ i 


(EGRID( 



( 1 ) 

!) ,NGRID(2) ,EGRID(1) ,NGRID( 1 
S),NGRID(3),EGRID(4),NGRID(4 



): 



AZ21.DIST21) 

AZ34.DIST34) 



COMPUTE APPROXIMATED COORDINATE OF TRAVERSE STATION POSITIONS 



READ NUMBER OF OBSERVED ANGLES 
READ(5,*) N 

READ NUMBER OF RESECTION POINT 
READ(5 ,*) NR1 



NSTA = N-l 
XGD = EGRIDC 2) 

YGD = NGRID(2) 

AZFIR = AZ21 

CALL DMSR (NUM360 ,NUM0 ,NUM0 .AN360R) 

K = 0 

DO 100 1=1, NSTA 

READ ANGl!E IN DEGREE, MINUTE, SECOND 
READ(5 , 1 0 5 j D.M.S 
Dd = DFLdAtfD 
MM = D FLOAT! M 
SS = D FLOAT (S 

CALL DMSR (DD,MM,SS,R) 



ANGl(I) = R 
“ ^(1) = 



STDAi 



DUM10 



READ REDUCE GRID DISTANCES 
READ (5 *) RANG 

DfST(I) = RANG 
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c 

c 

c 

c 

c 



READ(5,*) NZ1(I) 



IF (RANG. EQ. 0. 0) THEN 

READ APPROX. GRID COORDINATES OF RESECTION STATION 
READ( 5 ,*) NCO 
READ(5 *) ECO 



END IF 



CNGRID(I) = NCO 
CEGRID(I) = ECO 
GO TO 100 



K = K+l 

DIS(K) = DIST(I) 

STDD(K) = (DIS(K)*DUM11) + DUM12 

AZIMUT = AZFIR + ANGl(I) 

IF (AZIMUT. GT.AN360R) THEN 
AZIMUT = AZIMUT - AN360R 
END IF 



C 

cooooo 

c 



cooooo 

c 



COAZR(K) = AZIMUT 
FOR PRINTING 



CALL RDMS (AZIMUT,DUM1D a DUM1M,DUM1S) 
CFAZD(K) = 1 DINT(OUMID) 
CFAZM(K) = IDINT(DUMIM) 
CFAZS(K) = DUM1S 



CALL UTM (XGD.YGD, RANG, AZIMUT, DUMMY1.DUMMY2) 



DELX(K) = 
DELY(K) = 

SUMDX 

SUMDY 



C 

C 

C 



C0OOOO 

c 



CEGRID 

CNGRID 



( 1 ): 



DUMMY1-XGD 

DUMMY2-YGD 

SUMDX + DELX(K) 
SUMDY + DELY(K) 

DUMMY1 

DUMMY2 



C 

C 



CALL GRIDAZ (DUMMY1.DUMMY2, XGD.YGD, AZC.DUM1S) 
FOR PRINTING 

CALL RDMS (AZC,DUM1D,DUM1M,DUM1S) 

CBAZD(K) = IDINT(DUMID) 

CBAZM(K) = IDINT(DUMIM) 

CBAZS(K) = DUM1S 



AZFIR = AZC 
XCOX(K) = DUMMY 1 
YCOY(K) = DUMMY2 



C 

C 

C 

C 

100 CONTINUE 



XGD = DUMMY1 

YGD = DUMMY2 

TODIS = TODIS + RANGE 
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ooo 



c 

c 

c 



K10 = K^l 



9 

C 

C 

C 



c 



c 

c 

c 



c 



c 



c 

C@(3@(9@ 

C 

c 



c 

c 



C@(3(3(3(3 

c 



c 



101 

c 

c 



DO 9 I =1 , K10 

TTCOXm = XCOX(I) 
TTCOY(I) = YCOY(I) 
CONTINUE 



KK = K+l 

READ (5.105) D A M A S 
D t) = DFLdAtfD) 

MM = DFLOAT M) 

SS = DFLOAT(S) 

READ (5 *) RANG 

DiST(N) = RANG 

CALL DMSR (DD,MM,SS,R) 

ANG1(N) = R 
STDA(N) = DUM10 

AZLAS = AZFIR + ANG1(N) 

IF (AZLAS.GT. AN360R) THEN 
AZLAS = AZLAS - AN360R 
END IF 

DIFAZI = AZLAS - AZ34 

WPRED(1)= DIFAZI 

CORAZ = DIFAZI / DFLOAT(KK) 

FOR PRINTING 

CALL RDMS (DIFAZI ,DUM1D,DUM1M,DUM1S) 

EAZ1D = I D I NT( ABS (DUM1D) ) 
EAZ1M = IDINTC ABS (DUM1M) ) 
EAZ1S = ABS (DUM1S) 

CALL RDMS (C0RAZ,DUM1D,DUM1M,DUM1S) 

EAZ2D = I D I NT( ABS (DUM1D) ) 
EAZ2M = IDINTf ABS (DUM1M) ) 
EAZ2S = ABS (DUM1S) 



WPRED(2j = 



WPRED(3; 

DIFDIS 

ERR 



SUMDX - (EGRID(3) - EGRID(2) ) 
SUMDY - (NGRID(3) “ NGRID(2) ) 



= DSQRT ( 
= I DINT ( 



(WPRED(2))**2 + (WPRED(3))**2 ) 
TODIS / DIFDIS ) 



READ NUMBER OF MEASURED ANGLES AT EACH TRAVERSE STATION 



DO 101 J=1 , KK 
READ (5,*) N2 
nl)m(J) = 

CONTINUE 



N2 



DO 108 I = 1,NR1 
II = K10+I 
READ(5,*) DUMMY1 
READ(5!*j DUMMY2 
TTCOY(ll) = DUMMY1 
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108 

C 

C 

C 

C 

m 

m 

c 



c 

102 

C 

m 

c 

c 

C(3(3 

C 

C 

c 

c 

CPP 

C 

C 

CPP 

C 

C 

c 

C00 

C 

CP(3 

C 

C 



TTCOX(II) = DUMMY2 
CONTINUE ~ 

ASP = NSTA - 1 

KEPP APPROX. GRID COORDINATES OF EACH STATION IN 
ASNY AND ASEX 

DO 102 1=1, ASP 

ASNY(I) = CNGRID(I) 

ASEX(I) = CEGRID(I) 

CONTINUE 

NUMBER OF OBSERVED ANGLES PLUS OBSERVED DISTANCES 
NT = N+K 

NUMBER OF OBSERVED ANGLES 
NA = N 

NUMBER OF OBSERVED DISTANCES 
ND = K 

NUMBER OF TRAVERSE STATIONS NOT INCLUDED KNOWN STATIONS 
NP = ASP 

NUMBER OF UNKNOWN DX AND DY 
NP2 = K10*2 + NR1*2 

NI = NP2**2 + 3*NP2 

NUMBER OF TRAVERSE STATION (NOT INCLUDE KNOWN STATION) 
PLUSE NUMBER OF INTERSECTION POINT 

NC = K10+NR1 



C(3@(3(9(3 CALL SUBROUTINE TO ADJUST STATION POSITIONS 

C ‘ ‘ ‘ 

C 



C 

C 

C@@P@ 

c 



BY LEAST SQUARE METHOD OF OBSERVATION EQUATION 

CALL INDLSQ (NC.NR1,NZ1,NUM.KK,K.DIST.TTC0X JTCOY.NT.NA.ND.NP, 
NP2 , NI , AZ2l ,AZ34 A tfGfUD , EdRID A Aftal , STdA ,tflS , ST0D , 
PR I NTD PW , Wtf , STDO , FM , BM , BMT , &MTWM , NM , T$ , 

DELTA , NMI ,WKlO , ASNY , ASEX , VM , STDAD) 



WRITE(6,* 
WRITEi 6 * 
WRITE(6 * 



C 

C 



103 



PRINTING DETAILS 

| S|ANDARD DEVIATION OF UNIT WEIGHT =' , STDAD 
'STANDARD DEVIATION OF ADJ. ANGLES' 

DO 103 1=1. NA 

CALL RCTmS (VM( 1.1 ) , DUM1D , DUM1M , DUM1S) 

DUMMY 1 = 6UMiD*36O0. 0D0 + DUM1M*60. 0D0 + DUM1S 
WRITE(6 ,*) 'I =' , I , ADJ. STDA. =' , DUMMY 1 

CONTINUE ' 

WRITE(6,*)' STANDARD DEVIATION OF ADJ DISTANCES' 
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oooooooooooooo 



c 

c 



DO 104 1=1, ND 
K4 = NA+I 
WRIT€(6,*) 1 1 

WRITE(6 *) 'I = 1 , I , 1 ADJ . STDD. 



104 CONTINUE 



=' ,VM(1,K4) 



105 FORMAT (I3,1X ) I2,1X,I2) 

STOP 
END 



(3 (3 (3(3 (3 (3 (3 (3 (3 (3 (3(30(3(3(3(3(3 (3(3 (3 (3 (3 (3(3 (3 (3 (3 (3(3 (3 (3(3 (3(3(3 (3(9(3 


@@ 






@0 


@0 


SUBROUTINE LEAST S( 


1UARE 


@@ 


@@ 


ADJUSTMENT OF INDIf 


<ECT 


@@ 


@@ 


OBSERVATIONS 




@0 


@0 






@0 



@@(3(3 (3 (3 (3(3 (3 (30(3(3 (3 0(3 (3(9 (3 (3 (3(3(3 13(3(3(3(3 (3(3(3 (3 (3 (3 (3(3 (3 (30 



C 

c 



c 

c 

c 

c 



SUBROUTINE INDLSQ ( NCI , NR , NZZ , NUMM . KK, K , DI ST1 . TCOX . TCOY , NT , NA , 

* ND.tfPjJP2,tfl.AZ2l,Az34.NGRfD.EGGlD,Alta, 

* stOa.Oisi !stC)d,pr1ntd.?w,wm .Stdo.fm, 

* BM.BHT.BMtWM.NM.TM.DEtTA NMi.WKld, 

* ASfsIG ,A$ EG .VM; STAND) 

DOUBLE PRECISION DUMMY1 ,NGRID(4J ,EGRID(4) . DUM1D . DUM1M, 

* DUMlS,XNG(NAT,DfSlLND7.pj'NT > Nt).STD(}CNT), 

* ASNGf ftp} .ASEG(NP) ,FM(Nt) ,BM(flT,IW2) ,TESTl ! 

TEST2.BMtCNP2,NTi !BMTWM( nP2,NT) ,NM(NP2,NP2) , 
TM(NP2) .DELTA(NPz) ,NMI(NP 2 ,NP2) ’WKlO(Ni) , 
STDA(nA) , STDDf NDK6UMMY2 ,NUM0,A221 ,AZ34, 
VMC1, NT}! STAND, DI§T1(NA),TC0X(NC1), 

TCOY(NCl) 



I ,J ,K1, IER,PRINTD, PW. CHECK, NUMM(KK) ,K,KK, 

nf^, ^zz(na) , kon , i i i , nc!i ,nnd! 



c 

c@@ 

c 



_ INTEGER 

DATA K1/0/.TEST1/0. 00010000000D0/,NUM0/0. ODO/ 



NC4 = 4 

NP5 = 6 
N1 = 1 



20 

30 

C 

C@@ 

C 



SET WEIGHT MATRIX 

DO 30 1=1, NT 
DO 20 J=1 ,NT 

WM(I,j) = 0. ODO 
IF (i.EQ.JJ THEN 
WMil .JJ = 1. ODO 
END IF 
CONTINUE 
CONTINUE 

FOR STD. ANGLE 

IF (PW. NE. 1) THEN 
UO 40 1=1, NA 

DUMMY 1 = STDA(I) 
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40 



CONTINUE 



CALL DMSR (NUM0.NUM0.DUMMY1 .DUMMY2) 

STDO(I) = l.ODfl / ( (0SIN(DUMMY2) )**2) 



COO FOR STD. DISTANCE 
C 

DO 50 1=1, K 

J = NA+I , x 
DUMMY 1 = STDD(I) 

STDO(J)= l.ODO / (DUMMY1**2) 
50 CONTINUE 
C 

COO SET UNEQUAL WEIGHT 

DO 70 1=1, NT 
DO 60 J=1,NT 



60 

70 



M( 

END IF 
CONTINUE 
CONTINUE 

END IF 



IF (I.mjJ THEN 
WMil.JJ = STDO(I) 



C 

C 

< 

C' 

c 



COO PRINT WEIGHT MATRIX 
C 

IF (PRINTD.NE. 0) THEN 

WRITE(6,*) 'WEIGHT MATRIX' 

WRITEC6 *) ' v 

CALL USWFM ( ' R-C. ' ,NC4,WM,NT,NT,NT,NP5) 

END IF 

'99 CONTINUE 

K1 = Kl+1 
C 

IF (K1 .GT. 2 ) THEN 
GO TO 999 
END IF 
C 

CHECK = 0 
C 

C@@ CALL SUBROUTINE TO CALCULATE "F" MATRIX 

C 

C 

CALL CALFM (NUMM.KK.K.DIST1 .TCOX.TCOY ,AZ21 .AZ34.NA.ND, 

* NT.NtaRld.teRID.XNG.DfSl.N^ASNd.ASEd.FM) 

C 

c 

c 

c 

C@@ CALL SUBROUTINE TO CALCULATE "A" MATRIX (COEFFICEINT DX,DY) 

CALL CALAM ( NR, NZZ. NUMM.KK.K.DIST1 . TCOX.TCOY, NT .NA.ND.NP, 

* NP^,NG^ID,EdRld,ASNG,A^EG,BM) 

C 

C@@ TRANSPOSE OF "A" MATRIX 
C 

DO 120 1=1, NT 
DO 110 J=1 , NP2 

c BMT(J,d) = BM(I,J) 

110 CONTINUE 
120 CONTINUE 
C 

COO "A" TRANSPOSE * W 
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c 

CALL VMULFF (BMT,WM,NP2,NT,NT,NP2,NT,BMTWM,NP2, IER) 

"A" TRANPOSE * W * A 
C 

CALL VMULFF (BMTWM,BM,NP2,NT,NP2,NP2,NT,NM,NP2,IER) 

C3P "A" TRANSPOSE * W * F 
C 

CALL VMULFF ( BMTWM , FM , NP2 , NT , N1 , NP2 , NT , TM , NP2 , I ER) 

C 

c 

c 

CPP (INVERSE OF ("A" TRANSPOSE W * A) ) 

CALL LINV2F (NM.NP2 ,NP2 ,NMI ,N1 ,WK10, IER) 

C 

C(3(3 CALCULATE DELTA MATRIX (COORECTION VECTOR) 

C (INVERSE OF ("A" TRANSPOSE W * A) * "A" TRANSPOSE * W * F ) 
CALL VMULFF (NMI ,TM,NP2,NP2 ,N1 ,NP2,NP2 .DELTA ,NP2 , IER) 

C30 UPDATE APPROX. VALUE 

C 

C 

C KON = K-l+NR 



C 

C 



DO 130 1=1, NCI 

TCOX(I) = TCOX(I) + DELTA((I-1)*2 + 1) 
TCOY(I) = TCOY(I) + DELTA(( I-l)*2 + 2) 



130 CONTINUE 



III = (K-l)*2 
DO 131 1=1 ,NP 



I 



132 



F (NZZ(I) .EQ. 0) THEN 

ASEG(I) = ASEG(I) + DELTA( !I-1)*2 + 1 
ASNG(I) = ASNG(I) + DELTA( (I-l)*2 + 2) 

ELSE 

DO 132 J =1 , NR 

IF (NZZ(I) .EQ. J) THEN 

ASEG(I) = ASEG(I) + DELTA( II1+(J-1)*2 +1 
ASNGflj = ASNG(I) + DELTA( II1+(J-1)*2 +2 
GO TO 131 
END IF 
CONTINUE 
END IF 



131 CONTINUE 



IF (PRINTD .NE. 0) THEN 
WRITE! 6,1020) K1 

' CALL 3 3sWFM ( ' R-C. ' ,NC4 , FM , NT , NT , N 1 , NP 5 ) 
WRITE! 6 , 1000) 

* RITt CALL 4 8 sWFM ( 1 R-C. 1 ,NC4,BM,NT,NT,NP2,NP5) 
WRITE!6,1000) 

WRITE( MwFM ('R-C. ' ,NC4,BMT,NP2,NP2,NT,NP5) 
WRITE!6,1000) 

WRI "cALL 6 U5V/FM ('R-C. ' ,NC4, BMTWM, NP2.NP2, NT, NP5) 
WRITE!6,10001 

CAlL ? U$WFM ('R-C. ' ,NC4,NM,NP2,NP2,NP2,NP5) 



J 
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WRITE 6,1000 
WRIT£(6 1080} ,, 

CALL USWFM ('R-C. ' ,NC4,TM,NP2,NP2,N1,NP5) 
WRITE(6 ,1000) 

WRITEC 6 1090) 

CALL USWFM ( ' R-C. 1 ,NC4,NMI,NP2,NP2,NP2,NP5) 
WRITE(6 , 1000) 

WRITE(6 1100) 

CALL USWFM (R-C. 1 ,NC4, DELTA, NP2,NP2,N1,NP5) 
WRITE(6 , 1000) 

WRITE(6 ; 1110) 

DO 140 1=1 , NP 



140 



CONTINUE 
END IF 



WRITE(6*,*) I,' N = ' , ASNG(I) 

WRITE(6’*) f ' 



E = ' ,ASEG( I) 



C 
C 

C00 CHECK DELTA MATRIX 

C 

C 

DO 180 1=1, NP2 

TEST2 = DABS (DELTA(I)) 

IF ( TEST2 .GE. TEST1) THEN 
CHECK = 1 
END IF 

180 CONTINUE 
C 
C 

IF (CHECK .NE. 0) THEN 
GO TO 99 
END IF 
C 

COO CALCULATE "A" * X MATRIX 
C 

CALL VMULFF (BM,DELTA,NT,NP2,N1 ,NT,NP2,$TD0,NT, IER) 

C@0 CALCULATE "V" MATRIX 
C 

DO 190 1=1, NT 



STDO (f) = STDO(I) - FM(I) 
VM(1,I) = STDO(I) 

190 CONTINUE 



C 

C 

C 

C 



139 

149 



CORRECT OBSERVED ANGLES AND DISTANCES 

DO 139 I = 1,NA 

ANG(I) = ANG(I)+(STDO(I)) 

CONTINUE 
DO 149 1=1 . ND 
NND1 = NA+1 

DIS1(I) = D I S 1 ( I )+STDO(NNDl) 

CONTINUE 

IF (PRINTD . NE. 0 ) THEN 
WRITE(6 , 1000T 
WRITE(6. 1120) 

CALL USWFM ( 1 R-C. ' ,NC4,STD0,NT,NT,N1,NP5) 
WRITE(6 , 1000) ' 

WRITEC6 . 1130) 

CALL USWFM ( 1 R-C. 1 ,NC4,VM,N1,N1,NT,NP5) 

END IF 
C 
C 

C@@ CALCULATE "V" TRANSPOSE * W * V MATRIX 
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CALL VMULFF (WM,STD0,NT,NT.N1,NT.NT, FM.NT. IER] 
CALL VMtH-FF ( VM, FM ,Nl , Nt ,Ni ,Nl , Nt , StANO ,Ni , IER) 



C 

C 

C 

C 



IF (PRINTD . NE. 0) THEN 
WRITE(6 , 1000) 

WRITE(6 1140) 

CALL USWFM ( 1 R-C. 1 ,NC4 , FM,NT,NT ,N1 ,NP5) 
WRITE( 6 , 1150) 

WRITE(6 1160) STAND 
END IF 

CALL VMULFF (NMI , BMT.NP2 ,NP2 ,NT ,NP2 ,NP2 , BMTWM.NP2 , IER) 

STAND = DSQRT ( STAND / (DFL0AT(NT-NP2)) ) 

DO 240 1=1, NP2 
DO 230 J=1 ,NP2 

NMI (I,J) = STAND * (DSQRT(DABS (NMI(I.J)) ) ) 



END I 

'230 CONTINUE 
240 CONTINUE 



IF (I.EQ.J) THEN 
. lir ,TMir) = NMI(I,J) 



C 
C 

C@P CALCULATE "A"*( INVERSE "A"TRANSPOSE*W*A) * "A"TRANSPOSE 

CALL VMULFF (BM,BMTWM,NT,NP2,NT,NT ) NP2,WM,NT ) IER) 

IF (PRINTD . NE. 0) THEN 
wRITE(6, 1000 T 
WRITEC6 1170) 

CALL US&FM (‘R-C. 1 ,NC4,WM, NT.NT, NT.NP5) 



C 

C 



END IF 



DO 280 1=1, NT 
DO 270 J=1 .NT 

WM (I,J) = STAND * ( DSQRT ( DABS(WM( I , J) ) ) ) 



IF (I.EQ.J) THEN 
VMf 1 , 1 / = WM( I , J) 
END IF 



270 

280 

C 

999 
C 

C 

1000 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 



CONTINUE 

CONTINUE 



RETURN 



FORMA' 

FORMA- 

FORMA' 

FORMA- 

FORMA' 

FORMA' 

FORMA- 

FORMA' 

FORMA- 

FORMA' 

FORMA- 

FORMA' 

FORMA 

FORMA 

FORMA' 



' W1 sl 



1 

/// 



,5 
5X 
5X 
5X 
5X 
5X 
5X 
, 5X 
5X 
5X 
5X 
5X 



7 // 5X 
:///,5X 



1 A MATS T Y 1 /} 

'TRANSPOSE’OF A MATRIX' ,/) 

'A TRANSPOSE * W MATRIX'. 7) 

'A TRANSPOSE * W * A MATRIX',/) 
'A TRANSPOSE * W * F MATRIX 1 :/) 
'THE INVERSION OF A TRANSPOSE * 
'X MATRIX OR DELTA MATRIX',/) 
'UP-DATE THE APPROX. VALUES',/) 
'V MATRIX 1 ,/) , x 

'V TRANSPOSE MATRIX',/) 

'W * V MATRIX' ,/) 

'V TRANSPOSE * W * V MATRIX' ,/) 



W * A' ,/) 
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oooooooooooooo 



1170 mW-Mmto&E A TRANSPOSE *W*A)*A TRANSPOSE 1 ,/) 



END 



0 ( 5(8 000 

00@ SUBROUTINE FOR COMPUTING 000 
000 F MATRIX 000 
000 000 
000000000000000000000000000000000000000 



SUBROUTINE CALFM (K0UNT,KK,K^DIS } CQX,C0Y,AZ21 ,AZ34, 

* NNA , Nfto , MnY , GRV , GR)( ,OBANG , 060IS , 

* NOS ASY ASX F) 



C 

C 

C 

C 

C 

C 

C 

c 



c 

c 



c 

c 



c 

c 



DOUBLE PRECISION GRYjf 4) .GRX(4) .OBANGjfNNAl .OBDIS(NND) , 

N36u,n360R^DuMi , NUM u , CDl3( 100), TEMPI , 




I NTEGER 1 1 , 1 2 , K , KK , KKK , N , STN , KOUNT( KK) 

DATA N360/360. 0D0/ , NUMO/O. 0D0/ , N/0/ , STN/1/ 

CALL DMSR (N360,NUM0,NUM0,A360R) 

DO 500 1=1, NNA 
DISTAN = DIS(I) 

IF (STN . EQ. 1) THEN 
TEMPI = GRXm 



TEMP2 = GRY> 
AZB = AZ21 
TEMP3 = ASX' 
TEMP4 = ASY 



Hi 



CALL GRIDAZ (TEMPI ,TEMP2 ,TEMP3,TEMP4,AZF,DUM1) 

CDIS(I). = DUM1 
IF (AZF . GT,_ AZB1.THEN. 



ELSE 



F(I) = OBANG(T)-(AZF-AZB) 



FiT) = 0BANG(I)-(A360R+AZF-AZB) 



N = N+l 
GO TO 400 
END IF 



IF (STN .EQ. 2) THEN 
TEMPI = COa C STN— 



- VrSTN-1) 
IEMP2 = COY(STN-l) 



TEMPS = GRX(2) 

TEMP4 = GRY(2) 

CALL GRIDAZ (TEMPI ,TEMP2 ,TEMP3 ,TEMP4, AZB, DUM1) 
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ooooo 



IF 



C 

C 



C 

C 



c 

c 



c 

c 



c 

c 



c 

c 



( KK 

7EMP3 = GRX(3) 
TEMP 4 = GRY(3) 



EQ. 3) THEN 
' = GRX( ' 



END IF 

TEMP3 = ASX(I) 

TEMP4 = ASY(I) 

CALL GRIDAZ (TEMPI ,TEMP2,TEMP3 ,TEMP4,AZF,DUM1) 



CDIS(I) = DUM1 
" (AZF .GT. AZB) 

F(I) = OBANG(l )-(AZF-AZB) 



IF 
ELSE 



GT. AZB) THEN 



F(_I) = OBANG( I )-(A360R+AZF-AZB) 
END IF 



N = N+l 
GO TO 400 



END IF 



IF (STN .EQ. KK) THEN 
AZF = AZ34 



AZF = AZ3 
TEMPI = GRX 
TEMP2 = GRY 



IF (DISTAN . EQ. 0. ) THEN 
TEMP3 = AS X ( I ) 

TEMP4 = ASY(I) 

ELSE 

TEMP3 = COX(K-l) 

TEMP4 = COY(K-l) 

END IF 

CALL GRIDAZ (TEMPI ,TEMP2 ,TEMP3 ,TEMP4 , AZB, DIS1) 

IF (AZF .GT. AZB) THEN 

T(I) = OBANGCI)-(AZF-AZB) 

F^) = OBANG(I)-(A360R+AZF-AZB) 

N = N+l 
GO TO 400 
END IF 



IF (STN .EQ. K) THEN 
TEMPI = COXYK-i; 

TEMP2 = COY(K-l 
TEMP3 = COX(K-2 
TEMP4 = COY(K-2; 

CALL GRIDAZ (TEMP1,TEMP2,TEMP3,TEMP4,AZB,DUM1) 

IF (DISTAN .EQ. 0.) THEN 
TEMP3 = ASX(I) 

TEMP4 = ASY(I) 

ELSE 

TEMP3 = GRX(3) 

TEMP4 = GRY(3) 

END IF 

CALL GRIDAZ (TEMPI ,TEMP2 ,TEMP3 ,TEMP4,AZF,DUM1) 
CDIS(I) = DUM1 
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c 

c 



c 

c 



c 

c 



c 

c 



400 



C 

C 



c 

c 



IF (AZF .GT. AZB) THEN 

F(I) = OBANGCI)-(AZF-AZB) 
ELSE 



END I 

N = N+l 
GO TO 400 



F ( I ) = 0BANG( I )-(A360R+AZF-AZB) 



END 


IF 








IF 


(STN 




GT. 2) 


THEN 




TEM 


PI 


= COX( 


;STN-1) 




TEM 


P2 


= COY 


STN— 1 1 




TEM 


P3 


= COX 


'STN-2) 




TEM 


P4 


= COY( 


;stn-2) 




CALL i 


GRIDAZ 


(TEMPI 



TEMP3 = ASX(I) 

TEMP4 = ASY(I) 

CALL GRIDAZ (TEMPI ,TEMP2 ,TEMP3,TEMP4 , AZF, DUM1) 

CDIS(I) = DUM1 
IF (AZF . GT,_ AZB1 THEN. 



ELSE 



F(I) = OBANG(I )-(AZF-AZB) 



F(JJ) = OBANG(I)-(A360R+AZF-AZB) 
N = N + 1 



END IF 



IF ( (DISTAN .NE. 0. ) .AND. 
ADISISTN) = CDISm 
WRITE(6 , ) STN=‘ ,bTN , 



END IF 



KKK = KOUNT(STN) 

IF (N .LT. KKK) THEN 
STN = STN 
ELSE 

STN = STN+1 
N =0 
END IF 



(STN . LE. K) ) THEN 
CDI S= ' ,CDIS( I ) 



500 CONTINUE 



II = NNA 
DO 600 12=1, NND 
II = 11+1 

Ff II) = OBDI S( 1 2) - ADI S ( 1 2 ) 
600 CONTINUE 

RETURN 

END 
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555055555559055555555500505505559550500 



59 00 
59 SUBROUTINE FOR CALCULATING 99 
99 A MATRIX 99 
99 99 



555995095955959990999559999559999999999 



SUBROUTINE CALAM ( N R , NZ , COUNT , KK , K , DI SS , COXX , COYY , NT , NA , 
ND,NP,NP2,KllY,l(N)(,ANY,ANX,XM) 



DOUBLE PRECISION KNY(4) ,KNX(4)ANY(NP),ANX(NP_), DISTAN, 
AM^NT ,TjP2) iA DI^S( NA) , CdXX(K) , CQYY( K) , 



DUl ,DU2,Du3,DU4,DUS,DU6 



INTEGER 



1. 11,12,13,14 15 J1.J2, K, KK,K1,K2,K3, 

K4 ,N , NN , STfj ,T^M ,f EMl , T^m 2 , J j.jjl , jJ2 . 
NOM , COUNT(K^) , Tt , TT1 TEM3 , NU$ED , NZ( NA) , 
NR, L, 110 



DATA STN/1/ , N/l/ , NN/0/ , TT/0/ , J J/l/ 



NP1 = NP2-1 
K1 = NA-1 
K2 = K-i 
K3 = K+l 

110 = ( ( K-l ) *2) - 1 

(3 (3 (3 ( 3(3 (3 (3 (5 5(59 9 5 595559 559955 9 9 999 

5 9 

5 COEFFICEINT OF COORDINATES 9 
5 FOR ANGLE CONDITION 9 

9 5 

505505090955055555505505595559 



C 

C 

c 

c 

c 



231 

C 



C 

C 

232 



DO 800 I = 1 , NA 

DISTAN = DISS(I) 

DO 700 J1 = 1 ,NP1 ,2 
J2 = Jl+1 
AM(I,J1) = 0.0D0 
AMC I * J2) = 0. 0D0 

IF (STN . EQ. 1) THEN 



IF (DISTAN . EQ. 0) THEN 

DO 231 L=1 ,NR 

IF (NZ( I ) . EQ, L) THEN 
JJ = I10+2*L 
GO TO 232 
END IF 
CONTINUE 

ELSE 

JJ = 1 
i EM = JJ 
GO TO 232 
END IF 



IF (J1 . EQ. JJ) THEN 



ooo o o o o o oo o o ooo 



c 

c 



DU1 = ANY 
DU2 = ANX,_, . 

DU3 = (DUl**2) 

AM(I.Jl) = DU1/DU3 
AM(I,J2) = -DU2/DU3 




** 



2 ) 



N = N+l 
GO TO 700 
END IF 

GO TO 700 
END IF 

******************** 



IF (STN . EQ. 2) THEN 



IF 



(J1 

OU1 

DU2 

DU3 

DU4 

DU5 

DU6 



J) THEN 

(I) - COYY ( STN— 1 
(I) - COXXfSTN 
1**2)+(DU2**2) 






. EQ. JJ 
= ANY' 

= ANX 

= ^NY(2) _ - COYY(STH-l) 
= KNX(2) - COXX(STN-l) 
= (DU4**2)+(DU5**2) 



AM(I.Jl) = -(DU1/DU3)+(DU4/DU6) 
AM( I , J2) = (DU2/DU3)+C DU5/DU6) 



233 



IF (DISTAN .EQ. 0) THEN 
OO 233 L = 1 , NR 

IF ( NZ(f) . EQ. L) THEN 
JJ1 = I10+2*L 
GO TO 700 
END IF 
CONTINUE 

ELSE 

JJ1 = TEM + 2 
TEM1 = JJ1 
GO TO 700 
END IF 
GO TO 700 
END IF 

IF (J1 . EQ. JJ1) THEN 

AM(I,J1) = (DU1/DU3) 

AM( I ,J2) =-( DU2/DU3) 

N = N+l 
GO TO 700 
END IF 

GO TO 700 
END IF 

**************** 

IF ( (STN .GT. 2) .AND. (STN . LE. K2) ) 
K4 = STN - 2 



IF (J1 . EQ. JJ 
DU1 = ANY 
DU2 = ANX 
DU3 = (DU 



THEN 

I) - COYY(STN-l) 
I) - COXX(STN-l) 
t *2)+(DU2* > ‘"" 



DU4 = C0YY(K4) - COYY(Itn-I) 
DU5 = C0XX(K4) - COXX(STN-l) 



THEN 



78 



nno 



c 

c 

c 

c 

c 



234 

C 



c 

c 

c 

c 

c 

c 



c 

c 

c 



DU6 = (DU4**2)+(DU5**2) 



AM 

AM 



(I,J1) = -DU4/DU6 
(I,J2) = DU5/DU6 



JJ1 = TEM1 
GO TO 700 
END IF 

IF (J1 .EQ. JJ1) THEN 



M) : - (cul/a3) 



IF (DISTAN 



(DU2/DU3) - 



(DU4/DU6) 

(DU5/DU6) 



^ EQ. 0) THEN 

DO 234 .L = 1 , NR 



IF (NZ(I . EO. L) THEN 
JJ2 = I10+2*L 
GO TO 700 
END IF 



CONTINUE 
ELSE 

JJ2 = TEM1 + 2 
TEM2 = JJ2 
GO TO 700 
END IF 



GO TO 700 
END IF 

IF (J1 .EQ. JJ2) THEN 

AM(I.Jl) = DU1/DU3 
AM( I , J2) = -DU2/DU3 

N = N+l 
GO TO 700 
END IF 
GO TO 700 



END IF 



******************* 



IF (STN . EQ. K) THEN 
IF (J1 . EQ. Jj) THEN 
DUl = COYYfK-2' 

DU2 = COXX(K-2, 

DU3 = (DU1**2) + (DU2**2) 



DUl = DOYYCK-2) - COYY(K-l) 
“ ) - COYY(KtI) 



AM 

AM 



(I,J1) = -DU1/DU3 
(I,J2) = DU2/DU3 



JJ1 = TEM1 
GO TO 700 
END IF 

IF (J1 . EQ. JJ1) THEN 

IF (DISTAN . EQ. 0. ) THEN 
DU4 = ANY(I) - COYY( 

anx(i) - COXXf 

•- j4*^2) + (DU 



ELSE 



DU5 = 
DU6 = 



TEMP2 = J 



DU4 = 
DU5 = 
DU6 = 




K-l) 
K-l} 
5**2) 



- COYY( 

- COXXf 
2) + (DU 



K-l) 

K-l) 

5**2) 
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c 

c 



c 

c 
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C 

C 

C 

C. 

C 

c 

c 

c 

c 



700 



AM( I , Jl) = -(DU4/DU6) + (DU1/DU3) 
AM(I,J2) = (DU5/DU6) - (DU2/DU3) 

N = N + 1 

NUSED = Jl 

JJ2 = 0 
GO TO 700 
END IF 
NN = NN + 1 

AM( I , Jl) = -(DU4/DU6) + (DU1/DU3) 
AM(i;j2) = (DU5/DU6) - (DU2/DU3) 

DO 235 L = 1 .NR 

IF ( NZ( f) . EQ, L) THEN 
JJ2 = I10+2*L 
GO TO 700 
END IF 
CONTINUE 



GO TO 700 
END IF 

IF (Jl .EQ. JJ2) THEN 
AM( I , Jl) = DU4/DU6 
AM( I , J2; = -DU5/DU6 

N = N + 1 
GO TO 700 
END IF 
GO TO 700 



END IF 



****************** 



IF (STN . EQ, KK) THEN 
IF (Jl . EQ, JJ) THEN 

I)U1 = COYYfK-1) - KNY(3' 
DU2 = COXX(K-l) - KNX(3l 
DU3 = (DU1**2) + (DU2**^) 

AM( I , Jl) = -(DU1/DU3) 
AM(i;j2) = (DU2/DU3) 

N = N + 1 
GO TO 700 
END IF 



GO TO 700 
END IF 

CONTINUE 

NOM = COUNT(STN) 

IF (N . LE. NOM) THEN 
STN= STN 
ELSE 

STN = STN + 1 
JJ = TEM 
N = 1 
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= TEM] 



IF (STN 
~ JJ 

TEM1 = TEM2 
END IF 
END IF 



GT. 3^) THEN 



800 CONTINUE 



(? COEFF. OF COORDINATES 0 
@ FOR DISTANCE CONDITION @ 

(3(3(3(3(3(3(B(B(3(3(3(B(3(9(3(B(3(B(9(B(3(3(3P(9P 



C 

C 



c 

c 



c 

c 

c 



DO 3000 II = l.ND 

12 = NA + fl 

13 = II - 1 

14 = 13 + (13 - 1) 

15 = 13 + (13 + 1) 

DO 2000 J1 = 1.NP1.2 

J2 = J1 + 1 
AM(I2,J1) = 0. 0D0 
AM(I2 J2) = 0. 0D0 



IF (II 
IF 



.EQ. 1) 

VuVii 



THEN 
l) THEN 

1NYC2) - COYY(I) 
DU2 = KNXf2) - COXXfl) 
DU3 = DSQRT ( DU1**2 + 

AM( 12 , Jl) = — ( DU2/DU3 ) 
AM( 12 ’ v32 ) = -(DU1/DU3) 



DU2**2) 



END IF 
GO TO 2000 
END IF 



IF (II 
IF 



.EQ. ND) THEN 

(Jl . EQ. NUSED) THEN 

DU1 = COYY(K-I) - KNY(3) 

DU2 = C0XX(K-1] - KNX(3) 

DU3 = DSQRT ( DU1**2 + DU2**2) 
AM(I2,J1) = (DU2/DU3) 

AM(I2 J2) = (DU1/DU3) 



END IF 
GO TO 2000 
END IF 



IF 



(14 .EQ, Jl) THEN 
DUl = COYYCI3) - 



) - COYY(Il) 
) - COXX(Il) 
"001**2 + DU 



DU2 = COXXfl 3 
DU3 = DSQRT ( 

AM(I2,J1) = (DU2/DU3) 
AM( 12 ( J2 ) = (DU1/DU3) 

GO TO 2000 
END IF 

IF (15 .EQ, J1J THEN 
OU1 = COYYCI3) - COYY 
DU2 = COXX(.I3) - COXX 
DU3 = DSQRT ( DU1**2 



2**2 ) 



Hi) 

+ DU2* y 



2 ) 



i. 
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c 

c 



2000 



- AM( 12 , Jl) = 
AM( 12 | J2) = 

END IF 

CONTINUE 



-(DU2/DU3) 

-(DU1/DU3) 



C 

C 



3000 CONTINUE 
RETURN 

END 
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( 3 ^ 

@0 SUBROUTINE TO CALCULATE GRID 00 

00 AZIMUTH AND DISTANCE 00 

00 00 

( 3 ( 3(3 ( 3(3 ( 3(3 (3 (3 (3 (3 (3 (3 ( 3(3 (3 (3 (3 ( 3(3 ( 3(3 (3 (3 (3 (3 ( 3(3 ( 3(3 (3 (3 ( 3 ( 3(3 (3 ( 3 ( 3(3 

SUBROUTINE GRIDAZ (XI ,Y1 ,X2,Y2 , ANGLR , DIST12 ) 



REAL*8 X1.X2.Y1 , Y2.ANGLR.ANG0 .ANG90 . ANG180 .ANG270 
REAL*8 AN690ft,Afll8()R,AN2?0R,DiFX ,DlfY,DISTl2 

DATA ANG0/0. 0D0AANG90/90. 0D0/ .ANG180/180. 0D0/ 
DATA ANG270/270. 600/ 



CALL DMSR 
CALL DMSR 
CALL DMSR 

DIFX = X2-X1 
DIFY = Y2-Y1 



ANG90 ,ANG0 ,ANG0 .ANG90R ) 
ANG186 , ANGd , ANG6 , AN180R 
ANG270 ,ANG0 ,ANG0 ,AN270R 



IF((DIFX. EQ.ANGO). AND. (DIFY. EQ.ANGO)) THEN 
ANGLR = ANGO 

ELSE IFfDIFX. EQ.ANGO) THEN 
IF(DIFY.GT. ANGO) THEN 
ANGLR = ANGO 
END IF 

ANGLR = AN180R 



ELSE 



END IF 



IF(DIFX.GT.ANGO) THEN 

ANGLR = ANG90R-DATAN(DIFY/DIFX) 

ANGLR = AN270R-DATAN(DIFY/DIFX) 
END IF 



DIST12 = DSQRT(DIFX**2+DIFY**2) 

RETURN 

END 



(3(3(3000000000000000000000000000000000000 

00 SUBROUTINE FOR CHANGING THE 00 
00 ANGLE FROM DEGREE, MINUTE, AND 00 
00 SECONDS TO RADIAN 00 
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00 00 
00000000 ( 300000000000000000000000000 ( 3(300 

SUBROUTINE DMSR ( D1,M1,S1,R1 ) 

REAL*8 D1,M1,S1,R1,PI1 , N60 , N3600 , N180 

DATA N60/60. 0D0/ ,N3600/3600. 0D0/ ,N180/180. 0D0/ 

PI1 = 4. 0D0*DATAN(1. 0D0) 

R1 = ((D1+(M1/N60)+(S1/N3600))*PI1)/N180 

RETURN 

END 



000000000000000000000000000000000000000 
00 00 

00 SUBROUTINE FOR COMPUTING GRID 00 

00 COORDINATES FROM KNOWN AZIMUTH 00 

00 DISTANCE AND STARTING C00RDI- 00 

00 NATES 00 

00 00 

000000000000000000000000000000000000000 



SUBROUTINE UTM (XI ,Y1 ,DIS12,ANG12R,X2,Y2 ) 

REAL*8 X1,Y1,X2,Y2,DIS12 ,ANG12R,DIFX ,DIFY 

DIFX = DIS12*(DSIN(ANG12R)) 

DIFY = DIS12*(DC0S(ANG12R)) 

X2 = Xl+DIFX 
Y2 = Yl+DIFY 

RETURN 

END 



C 

c 

c 



c 



000000000000000000000000000000000000000 
00 00 

00 SUBROUTINE FOR CHANGING THE 00 

00 ANGLE FROM RADIANS TO DEGREE, 00 

00 MINUTES, AND SECONDS 00 

00 00 

000000000000000000000000000000000000000 

SUBROUTINE RDMS ( RR1,D2,M2,S2 ) 

REAL*8 RR1,D2,M2,S2,PI1,N60,N180,TDEG,DIF,TMIN 

DATA N60/60. 0D0/ ,N180/180. 0D0/ 

PI1 = 4. 0D0*DATAN( 1. 0D0) 

TDEG= (N180*RR1)/Pll 
D2 = DFLOAT(IDINT(TDEG)) 

DIF = TDEG-D2 
TMIN= DIF*N60 
M2 = DFLOAT(IDINT(TMIN)) 

DIF = TMIN-M2 
S2 = DIF*N60 

RETURN 

END 
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